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By 
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Ravinder Bhati^ and Cai Lin 2 
ABSTRACT 

Parallel computing studies are presented for a variety of structural 
analysis problems. Included are the substructure planar analysis of 
rectangular panels with and without a hole, the static analysis of space 
mast, using NICE/SPAR and FORCE, and substructure analysis of plane 
rigid-jointed frames using FORCE. The computations are carried out on the 
Flex/32 Multicomputer using one to eighteen processors. The NICE/SPAR 
runstream samples are documented for the panel problem. For the 
substructure analysis of plane frames, a computer program is developed to 
demonstrate the effectiveness of a substructuring technique when FORCE is 
enforced. On-going research activities for an elasto-plastic stability 
analysis problem using FORCE, and stability analysis of the focus problem 
using NICE/ SPAR are briefly summarized. Speedup curves for the panel, the 
mast, and the frame problems provide a basic understanding of the 
effectiveness of parallel computing procedures utilized or developed, 
within the domain of the parameters considered. Although the speedup 
curves obtained exhibit various levels of computational efficiency, they 
clearly demonstrate the excellent promise which parallel computing holds 
for the structural analysis problems. 
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SUBSTRUCTURE ANALYSIS OF RECTANGULAR PANEL WITH HOLE 


USING NICE/ SPAR AND FORCE 


1.1 Problem Description 

A two-dimensional elastic stress analysis of a rectangular panel with 
a central hole is conducted using NICE/SPAR while invoking concurrent 
processing with the use of FORCE. References 1 through 5 represent the 
relevant documents on NICE/SPAR and FORCE. Figure 1 shows the rectangular 
panel dimensions and the applied uniform compressive loading. The Young's 
modulus and the Poisson's ratio of the panel material are, respectively 

10,000 ksi, and 0.3, which represent a typical aluminum alloy. The panel 

thickness is taken as 0.1 in. Initially, the panel is analyzed in the 
absence of a central hole using substructuring. The analysis of panel 
including the hole is then conducted using two substructures representing 
the regions around and away from the hole. 

1.2 Finite Element Discretization and Substructuring 

Figures 2 through 4 show the various types of discretizations 

investigated using E41 quadrilateral elements of NICE/SPAR. Figure 2(a) 

shows the rectangular panel divided into four elements as well as the 
boundary simulations. Figure 2(b) shows substructures 1 and 2 with 
revised element and node numbering at the substructure level. Similarly 
Figure 3 shows an eight-element discretization and the corresponding 
substructures. Figure 4 shows the finite element discretization of the 
rectangular panel with hole. Figure 4(a) shows the node numbering scheme 
commencing from the inner region around the hole for nodes 1 through 24, 
and continues at the bottom left corner of the panel with nodes 25 




through 78. Figure 4(b) shows the element numbering scheme commencing 
from the inner region around the hole for elements 1 through 16, and for 
the remainder of the panel for elements 1 through 44 starting at the 

bottom left corner. The common element numbering from 1 to 16 can be used 
as long as the outer and the inner regions of the panel are modelled as 
separate substructures as shown in Figures 5(a), and 5(b), respectively. 

1.3 Solution Procedure and Parallel Computations 

The solution procedure for the substructuring method using NXCE/SPAR 
is given in this section. The relevant SPAR processors and CLAMP 
procedures. The steps of the procedure are as follows: 

1. The AUS subprocessor TABLE is created to input the number of 

substructures, the number of nodes per element, and the number of 
elements in each substructure. 

2. Interface nodes are entered using subprocessor TABLE. These are 

stored in the NICE/SPAR library in the form of data sets. 

3. Steps 1 and 2 are stored in a procedure called SUB_TABLE. 

4. An external file PRESERVE is introduced into the main runstream. 
This file consists of CLAMP procedures ENDJSUB and GLOB_DAT. 

5. The Individual substructure data such as the number of nodes in the 

substructure, material properties, joint locations, boundary 
conditions, and element connectivity are entered. At the end of each 
substructure data, procedure ENDJSUB is called to store the data set 
sequence number. 

6. Procedure GLOB DAT is called for generating the data sets required 

for the assembled structure, from the data obtained from each 
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substructure 


7. Processor K is executed to obtain the final stiffness matrix. 

8. subroutine getrow is used to extract the K matrix from the NICE/SPAR 
data library, and stored by rows in the upper triangular part of the 
matrix. 

9. The FORCE subroutine SGEFA is called to perform the decomposition of 

the square matrix using gaussian elimination with partial pivoting. 

10. Subroutine SGESL is called to solve the system of equations for the 

static displacements by back substitution. 

11. The FORCE subroutines in steps 9 to 10 are run on several processors 
for speedup study. 

1.4 Numerical Results and Discussion 

In this section as well as sections 2 and 4 of this report, the 

following definitions of speedup and efficiency are used. Speedup is 

defined as the execution time on a uniprocessor divided by execution time 
on n processors. Parallel efficiency is defined as the speedup divided by 
the number of processors times 100. Table 1 presents the speedup and 
efficiencies for the four element rectangular panel shown in Figure 2. 
The corresponding speedup versus the number of processors relationship is 
shown in Figure 6, involving a 27 x 27 stiffness matrix. The maximum 
speedup factor is 3.7 for 8 processors. The efficiency reduces 

continuously with an increase in the number of processors. 

Table 2 presents the speedup and efficiencies for the eight-element 
rectangular panel shown in Figure 3. The corresponding speedup versus 
number of processors relationship is shown in Figure 7, involving a 45 x 
45 stiffness matrix. The maximum speedup factor is 5.14 for 8 and 12 
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processors. The efficiency reduces continuously with an increase in the 
number of processors, though at a lower rate than that for the four 
element panel. 

Table 3 presents the speedup and efficiencies for the rectangular 
panel with a central hole shown in Figure 4. The corresponding speedup 
versus the number of processors relationship is in Figure 8, involving a 
234 x 234 stiffness matrix. The maximum speedup factor is 13.2 for 18 
processors. Here the speedup increases continuously with an increase in 
the number of processors while the efficiency decreases gradually* 
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2. CONCURRENT STATIC ANALYSIS OF MAST USING FORCE 


2. 1 The Mast Problem 

Figure 9 shows the mast problem. The mast is a 60-meter high truss 
divided into 55 platforms with three nodes at each platform level. Each 
node has three degrees of freedom, that is, displacements in x, y, and z 
directions. The mast consists of horizontal batten members, inclined 
diagonal members, and vertical longeron members. The bottom platform with 
first three nodes is constrained against displacements. A concentrated 
load of 100 tons is applied at the upper most node in the global x 
direction. The datasets were stored in the NICE/SPAR library named 
MASTM3.L01* The problem is to determine the static nodal deflections 
using NICE/SPAR and FORCE. 


2.2 Results and Discussion 

The stiffness matrix generated by NICE/SPAR is of the order 495 X 
495, including the support nodes. The linear system of equations is 
solved using FORCE. Table 4 presents the speedup and efficiencies when 1, 
2, 4, 8, 12, 16, and 18, processors are used to solve the system of 
equations. The highest speedup is obtained when 18 processors are used, 
and is 15.71. The computational efficiency decreases slightly with an 
increase in number of processors. The results are shown by the dashed 
curve in Figure 10, in which the theoretical ideal relationship is also 
shown for a direct comparison. In comparison to the results presented in 
Section 1 for the rectangular panel problem of the order 234 X 234, the 
speedups for the mast show better promise due to an increase of the 
stiffness matrix size to 495 X 495. 
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3. ELASTO-PLASTIC STABILITY ANALYSIS USING FORCE 


3. 1 Biaxially Restrained Imperfect Column 

In a study reported in Reference 6 by Darbhamulla, Razzaq, and 
Storaasli, a concurrent elasto-plastic solution for the stability analysis 
of biaxially restrained imperfect column was presented which employed the 
Finite Element Machine. Appendix C presents the listing of the 
corresponding computer program in PASCAL. This program is being converted 
to CONCURRENT FORTRAN so that the efficiency of the algorithm can be 
studied using FORCE developed by Jordan and Norbert (Reference 5) and 

implemented on FLEX/32. As described in Reference 6, the problem 
involves solving three coupled nonlinear differential equations whose 
coefficients vary with the applied loads and spatial coordinates. In that 
sense, the concepts developed are of general use in various types of 
materially nonlinear structural problems. A brief outline of the algorithm 
utilizing FORCE is given in the following section. 

3.2 Concurrent Elasto-Plastic Algorithm Using FORCE 

The computer program given in Appendix C is being converted to 
CONCURRENT FORTRAN based on the following algorithm: 

1. Read initial data and compute required cross-sectional properties. 

2. Assemble a global stiffness matrix and a force vector concurrently 
uing FORCE . 

3. Compute initial displacement vector using the FORCE simultaneous 
equation solver. 

4. Synchronize all processors. 

5. Compute elasto-plastic cross-sectional properties at various 
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locations along the length concurrently using the procedure 
presented in Reference 8. 

6. Reassemble the global stiffness matrix and the new force vector 
concurrently, and update the displacement vector using, again, the 
FORCE simultaneous equation solver. 

7. Check convergence and proceed to Step 8 if convergence is 

achieved. Otherwise, go to Step 4. 

8. Increment the external load and go to Step 4 if the member 
collapse has not occurred. 

The program is being implemented on the FLEX/32 MMOS computer. 
3.3 Publication 

Reference 8 is scheduled to appear in Engineering with Computers, An 
International Journal for Computer-aided Mechanical and Structural 
Engineering this year. The proof-copy of this paper is given in 

Appendix D. 
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4. SUBSTRUCTURE ANALYSIS OF PLANE FRAMES USING FORCE 


4.1 Problem Definition and Goals 

The problem is to develop a concurrent algorithm for the linear 
load-deflection analysis of rigid-jointed plane frames using 
substructuring. Figure 11 (a) shows an example of a three story frame 

ABCDEFGH with a span L and story height h. Figure 11 (b) shows one way of 
dividing the structure into two substructures, CDEF and ABCFGH. The main 
goal of this study is to evaluate the computational efficiency of the 
concurrent algorithm developed while employing FORCE when substructures 
of the type shown in Figure 11 are adopted. For the present, a special 
purpose matrix stiffness substructure analysis program has been developed 
and implemented on FLEX/32 parallel computer for a direct investigation of 
the effectiveness of FORCE, as the main objective. The concurrent 
substructuring technique being considered may be generalized for arbitrary 
frames for implementation on NICE/SPAR later. 

4.2 Substructuring Technique 

The substructuring technique adopted for the analysis of the frame 
considered here is given in Reference 9 and summarized in Appendix E of 
this report. A number of equations presented in Appendix E are used in a 
concurrent algorithm outlined in Section 4.3. The following is a brief 
summary of the applicable equations. The various terms are defined in 
Appendix E. 

The static joint equilibrium matrix equation for the k-th 
substructure is given by Equation E.l as follows: 

[S]{A} = {P} + {R} (!) 
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Equation E.4 for the k-th substructure can be written as: 


' s uAVk ■ !P A- 


( 2 ) 


for k = 1, 2, 3,..., n. The reactions at the restrained nodes can be found 
using Equation E.5: 


{R r } k = fS ru ] k^u } k’ 


(3) 


since {P f } - {0} for the substructures shown in Figure 11 (b) owing to 

the fact that the nodes C, F, A and R, are completely restrained. 

The substructure stiffness matrix given by Equation E.12 can be 
expressed as: 


tS s*k " ^rr^k'^nA^mA ^uA 


(4) 


If a concurrent inversion subroutine is unavailable, the matrix 
- 1 . 


product [S ]. [S ]. can be obtained as follows: 
r uu k ur k 

[A - [S ff [S ], 
ur k uu k l ur J k 

which may be expressed as: 

[S mA [A uA " ^ S ur^k 

where: 


^ ur ^ k ~ t* A ul* * A u2*’“* X urA 
fS ur ] k = [{S ul } {S u2 } ** ,{S ur }] k 


(5) 

( 6 ) 


(7) 

( 8 ) 


and {} are column vectors of the order u x 1. Thus, Equation 6 can be 
written as a series of matrix equations as follows: 

lS mA* A ul*k " {S ul } k 
fS mA {A u2 } k " 


u2 k 


(9) 

[S mA*\A " {S ur } k 

Solving the system of Equation 9 for {* ul > k » {A u2 } k »*« { A ur } k 

by means of a simultaneous equations solver, the matrix [ A 
Equation 6 can be formed. This clearly provides the matrix product 


Q 



indicated in Eq. 4 without inverting the matrix [ s uu l k * Equation 4 
can now be written as: 


' S A - < S rA-' S ruV\r’k 


( 10 ) 


The static equilibrium equation for the nodes linking the various 
substructures is given by Equation E.13 as follows: 


IS a ]fc a ( - (P a ) + [R a l, 


( 11 ) 


in which the various terms of [S ] and {P } are obtained as explained 

cl 3 . 

in Section E.3 of Appendix E. 

The true deflection vector, {a for the nodes linking the 

1U 

various substructures is found using Equation E.16 in the following form: 


< P rd> ' < S lll<\d> 


( 12 ) 


The deflections of the unrestrained nodes for the k-th substructure is 
found using Equation E.19. 


<\n>k = «V-t S ur^rd»k 


(13) 


The procedure described in this section can be readily applied to any 
type of framed structures behaving in a linear elastic manner. A 
concurrent algorithm utilizing this technique is given in the following 
section. 


4.3 Concurrent Algorithm 

Based on the substructuring technique outlined in Section 3.2, the 
following concurrent algorithm has been developed on FLEX/32 parallel 
computer. 

1. Formulate the stiffness matrix [S] given in Equation 1 for the 
substructures k *■ 1, 2, 3,..., m, on n processors concurrently. 

2. Extract [S ], appearing in Equation 2, from the [S] matrices, 

concurrently, for all k values. 
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3. Solve Equation 2 for {a u }^ with k = 1, to determine the 

deflections of Substructure 1 by invoking FORCE simultaneous 
equations solver SOLVE, concurrently using n processors. Similarly, 
solve Equation 2 for {A u > k with k = 2, 3,..., m, corresponding 
to Substructures 2, 3,..., m, using FORCE. 

4. Extract [S ], appearing in Equation 3, from the [S] matrices, 

ru k 

concurrently for all k values. 

5. Compute the vectors {Rj.}^ with k ■ 1 to determine the reactions 

in Substructure 1 by invoking FORCE matrix multiplication subroutine 
called MATMP, concurrently. Similarly, determine {R^^ for 

k =• 2, 3,..., m, corresponding to Substructures 2, 3,..., m, using 

FORCE. 

6. For k =* 1, determine {A with v = 1 as represented in 

Equations 9 using SOLVE, concurrently on n processors. Repeat the 
process of finding { A for v - 2, 3,..., r, respectively, to 

generate the matrix ^ ur^k as indicated in Equation 7. Obtain 
[S ], using Equation 10 concurrently on n processors. 

7. Repeat Step 6 for k - 2, 3,..., m. 

8. Assemble the global matrix [S ] and the vector {P } appearing in 

a a 

Equation 11, on a single processor using the [S^^ matrices 

generated in Steps 6 and 7; {R f } k generated in Step 5, while 

employing Equations E.14 and E.16 given in Appendix E. 

9. From [S ] and {P } obtained in Step 8, formulate Equation 12. 

3 3 . 

10. Solve Equation 12 concurrently, using SOLVE, on n processors. 

11. For k = 1, determine the true deflection {A } of the 

un iv 

unrestrained nodes, using Equation 13, concurrently on n processors. 

12. Repeat Step 11 for k ■ 2, 3 
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The algorithm outlined above is coded in CONCURRENT FORTRAN for the 
problem shown in Figure 11. A listing of the program is given in Appendix 
F with a sample input /output. The computational efficiency of the 
algorithm is explained in the following section. 

4.4 Numerical Results for 3-Story Frame 

With three degrees-of-freedom per joint (vertical and horizontal 
displacements, and a rotation), the frame stiffness matrix is of the order 
18 x 18. Table 5 presents the computational time t, speedup factor s, and 
efficiency for 1 to 5 parallel processors for a concurrent analysis of the 
3-story frame. The time t is measured in rtick units. One rtick equals 
l/50th of a second. For solving the 18 simultaneous equations 
concurrently, Subroutine SOLVE developed by Jordan and his research 
associates from the University of Colorado, was utilized. A complete 
listing of the subroutine appears as a part of the computer program 
presented in Appendix F. The speedup factor versus the number of 
processors data given in Table 5 is used to plot the 'actual' curve in 
Figure 12 in which the theoritical ideal relationship is also given for a 
comparison. The maximum speedup is 2.60 and is obtained when a total of 4 
processors are used corresponding to a computational efficiency of 65%. 
The use of 2 processors, however, results in an efficiency of 81%. 

Table 6 presents the computational times T^ through T^, defined 
as follows: 

T^ * time taken to assemble the stiffness matrices [S] given in 

Equation 1, for Substructures 1 and 2 shown in Figure 11, 

T 2 = time taken for solving Equations 2, 3 and 4 for Substructure 1 

with each of the matrices S , S and S of order 6x6 

uu ’ ru s 


12 


time taken for solving Equations 2, 3 and 4 for Substructure 


2, with the matrices S , S and S of the order 6x6, 

* uu * ru s 

12 x 6 and 12 x 12 respectively, 

=» time taken for solving Equation 12 for the complete system; 
the speedup factor, s, based on the total time T given by: 

T t * Tj + X 2 + T 3 + T 4 ; 

and efficiency based on s; for 1 to 5 parallel processors for a 

concurrent substructure analysis of the 3-story frame. The FORCE 
subroutine SOLVE was used concurrently for the two substructures for 
obtaining solutions to a pair of six simultaneous equations. It is noted 
here that there are two levels of computational parallelization involved 


in this solution scheme: 


Level 1: parallel generation of the two substructure stiffness matrices 
on two different processors, and 

Level 2: parallel solution of two sets of six simultaneous equations 

using SOLVE on up to five processors. 

The parallel processors work on one set of equations at a time. 

The value of T^ ranges from 3 to 5 rticks. The T^ value is in the 

range from 3 t 4 rticks, whereas T^ is in the range from 5 to 10 

rticks. One may expect to obtain nearly identical values for T 2 and 

T^ since they represent the computational times for handling the two 

substructures with equal number of free nodes. That this is not indeed 

the case can be understood by realizing that the orders of S and S 

ru s 

matrices for Substructures 1 and 2 are not equal, as described above in 


the definitions 

for 

T 2 and 

T 3* 

The 

value of T. is halved when 
4 

2 to 5 

processors are 

used 

instead 

of 

just 

one. The total time T 

varies 

between 12 and 

21 

rticks , 

the 

latter being the time taken 

by one 
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processor. For 2, 3 and 4 processors, T value is nearly constant. For 

5 processors it is greater than for 2, 3 or 4 processors. The maximum 
speedup is 1.75 and is obtained when a total of 3 or 4 processors are 
used. However, the efficiency with 3 processors is greater than that with 
4 processors. The maximum efficiency is obtained with 2 processors and is 
80.5%, of course excluding the case of one processor. 

The reason for the relatively low speedups obtained is that the 
matrix sizes for the problem considered are relatively small. As the 
number of processors is increased, the overhead time (for example, while 
executing the Barrier, End Barrier and Join statements) becomes 

significant as compared to the time consumed by arithmetical operations 
required for the simultaneous equations. 
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5. BUCKLING ANALYSIS OF FOCUS PROBLEM USING NICE/SPAR 


The buckling solution to the CSM Focus Problem 1 outlined in 

Reference 2 is being run for various discretizations. The finite element 
meshes are generated using CSM1 with NICE/SPAR E43 finite elements. 
Presently, all computations are being carried out on a single FLEX/32 
processor. The implementation of concurrent buckling and postbuckling 
solution procedures will be carried out once rigorous study of these 
problems is completed on a single processor. 
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Table 1. Speedup and efficiencies for four-element rectangular panel 


Number 

of 

Processors 

Computational 

time 

(secs) 

Speedup 

Efficiency 

1 

0.96 

- 

- 

2 

0.54 

1.78 

89.0 

4 

0.28 

3.43 

85.7 

8 

0.26 

3.70 

46.2 

12 

0.30 

3.20 

26.7 

16 

0.42 

2.29 

14.3 

18 

0.48 

2.00 

12.0 
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Table 2. Speedup and efficiencies for eight— element rectangular panel 


Number 

of 

Processors 

Comput a t io nal 
time 
(secs) 

Speedup 

Efficiency 

i 

3.60 

- 

- 

2 

1.88 

1.91 

95.7 

4 

1.06 

2.25 

56.2 

8 

0.70 

5.14 

64.3 

12 

0.70 

5.14 

42.8 

16 

0.84 

4.30 

27.0 

18 

0.92 

3.91 

21.7 
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Table 3. Speedup and efficiencies for rectangular panel with central hole 


Number 

of 

Processors 

Computational 

time 

(secs) 

Speedup 

Efficiency 

1 

474.62 


- 

2 

248.74 

1.90 

95.4 

4 

128.82 

3.69 

92.1 

8 

66.80 

7.10 

88.8 

12 

47.50 

9.99 

83.3 

16 

38.56 

12.31 

77.0 

18 

35.96 

13.20 

73.3 
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Table 4. Speedup and efficiencies for the mast problem 


Number 

of 

Processors 

Computational 

time 

(secs) 

Speedup 

Efficiency 

1 

4530.02 

- 


2 

2388.80 

1.90 

94.8 

4 

1228.10 

3.69 

92.3 

8 

610.00 

7.43 

93.0 

12 

417.38 

10.85 

90.5 

16 

320.60 

14.13 

88.3 

18 

288.30 

15.71 

87.3 
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Table 6. Computational time measurements, speedup factors, and 
effciencies for concurrent substructure analysis of 
3-story frame 


Number o f 

Execution time 

t (rticks) 

Total 

Time 

Speedup 

Factor 

Efficiency 

*? 

% 

Processors 

T 1 

T 

2 

T, 

mm 

1 

5 

4 

10 

2 

21 

1.00 

100.0 

2 

3 

3 

6 

i 

13 

1.62 

80.5 

3 

3 

3 

5 

i 

12 

1.75 

58.3 

4 

3 

3 

5 1 

i 

12 

1.75 

43.7 

5 

4 

4 

6 

i 

15 

1.40 

28.0 
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30.0 in 


Figure 1. Rectangular panel with hole and uniform compression 
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Figure 3. Substructuring of panel with eight elements 
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(a) Node numbering scheme 



(b) Element numbering scheme 


Figure 4. Finite element discretization of panel with hole 
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(a) Substructure 1 




(b) Substructure 2 


Figure 5. Substructuring of panel with hole and sixty elements 
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Speedup factor versus number of processors for 
four-element rectangular panel 
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Number of Processors 


Speedup factor versus number of processors for 
eight -element rectangular panel 
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Figure 11. Framed structure with two substructures 
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Figure 12. Speedup factor versus the number of processors for 
3-story frame using FORCE Subroutine SOLVE 
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Figure 13. Speedup factor versus the number of processors for 
3-story frame with two substructures using FORCE 
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APPENDIX A 


RELEVANT SPAR PROCESSORS AND CLAMP PROCEDURES 
A. 1 Relevant SPAR Processors 

A number of SPAR processors are described in Reference 3 some of 
which are used in the solution procedure presented in section 1.3 of this 
report. A brief description of the relevant SPAR processors is given in 
this section. The corresponding page numbers from volume 1 of reference 3 
are indicated in the parentheses following the SPAR processor names. 

1. Processor TAB (3.1-1) 

The function of this processor is to define the finite element model 
of the structure. The following subprocessors of TAB are utilized: 

JLOC defines the joint locations in a structure; 

MATC defines the material constants of the structure; 

JREF defines the orientation of reference frame associated with the 
joints ; 

CON is used in defining constraints and applied loading; and 
SA generates a table of shell (or panel) section properties to 

which reference is made during the definition of SPAR element 
E41. 

All the subprocessors generate their respective data sets which are 
stored in a data library. 

2. Processor ELD (3.2-1) 

This processor produces datasets containing element definitions such 
as the connecting joints, the type of SPAR element used, and the integers 
pointing to applicable lines in the table of section properties. 
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3. Processor AUS (5.1-1) 

This processor is comprised of an array of subprocessors which are 
data set constructors. The subprocessors are summarized in tab 5-1 (page 
5. 1-2, reference3) according to their functional categories. The 

functional categories include matrix arithmetic and modification of data 
tables. The following subprocessors of AUS are utilized: 

TABLE (2.5.1) creates the data sets which details the number of 
interface nodes, the number of substructures, and the number of 
elements in a substructure; and 

SYSVEC (5.1.3) is used primarily to represent either the joint 
displacements and rotations, or the applied forces and moments. 

4. Processor E (3.3-1) 

This processor generates in a skeletal form an element information 
packet for each element in the structure, and supplies general information 
such as the connecting joint numbers, material constants, and section 
properties data previously identified by processor ELD. 

5. Processor EKS (3.4-1) 

This processor reads data sets created by processor E and generates 
the element stiffness matrices. 

6. Processor TOPO (4.1-1) 

This processor analyzes element interconnection topology and creates 
data sets (KMAP and AMAP), Reference 4) for the assembly of the system 
stiffness matrix in SPAR standard sparse-matrix format. 

7. Processor K (4.2-1) 
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This processor assembles unconstrained system stiffness matrix in a 
SPAR standard sparse-matrix format. The input for K comes from the 
processors EKS, TOPO and partly from TAB. The output is a data set named K 
SPAR containing the unconstrained system stiffness matrix 

8. Processor INV (4.5-1) 

This processor factors assembled system matrix in the SPAR standard 
sparse-matrix format. It requires the unconstrained system matrix, the 
constraint definitions, and the topological information in AMAP. 

9. Processor SSOL (6.3-1) 

This processor computes the displacements and reactions due to applied 
loading at the joints. 

10. Processor GSF (7.1-1) 

This processor generates data sets containing element stresses. The 
input data sets are structural deformations from SSOL. The output data 
sets are named STRS. 

11. Processor PSF (7.2-1) 

This processor prints the element stresses and stress resultants from 
the data sets generated by GSF. 

12. Processor VPRT (5.3-1) 

This processor is used to edit and display information in SYSVEC 
format, such as the static displacements and reactions. It is also called 
the Vector Printer. 
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A. 2 CLAMP Procedures 


The substructure analysis using NICE/SPAR requires the use of MACRO 
subprocessors and CLAMP procedures described below. 

1. MACRO Subprocessor 

This is a subprocessor within the processor AUS, and is used to define 
a macrosymbol. The macrosymbols are used to establish variable like items 
which can perform arithmetic manipulations, access built-in functions, and 
rename directives. 

2. PROCEDURES 

The substructure CLAMP procedures END SUB and GLOB_DAT are used to 

minimize the user input data. These two procedures are stored in a data 
file named PRESERVE in the runstream presented in Appendix B. The 
procedures are described briefly as follows: 

Procedure END SUB scans through the table of contents and stores the 
sequence number of data sets. Later, the sequence numbers are used for 
global topology generation. The data sets are then renamed to make them 
distinct in order to avoid disability of the existing data set. 

Procedure GLOB DAT consists of a number of intermediate procedures 
such as Inter, Glob Jloc, Glob Mate. Procedure Inter searches the 
interface table for interface joints and separates it. procedure Glob_Jloc 
translates the substructure joint numbers to global joint numbers, and 
sets up global joint location data set. Procedure Glob_Matc creates the 
data set of material constants for the assembled structure. 
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APPENDIX B 


NICE/ SPAR RUNSTREAM SAMPLES 

The runstreams for the eight-element rectangular panel with hole are 
presented in this appendix including the outputs. The attached computer 
listings include: 

1. Runstream for eight-element rectangular panel. 

2. Output for eight-element rectangular panel. 

3. Runstream for rectangular panel with hole. 

4. Output for rectangular panel with hole. 

A listing of the FORCE program for concurrent solution of simultaneous 
equations developed by Jordan and his research associates, from the 
University of Colorado, is also contained in item 3 mentioned above. 
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RUNSTREAM FOR EIGHT-ELEMENT RECTANGULAR PANEL 


rmprocl 
rm SUBJABLE 
rm SS1.L01 

nicespar <<\endinput 
*set echo=off 
*open 1 ssl.101 /new 


DUM4Y START CARD 
[xqt tab 
start 1000 
1 

•fcnrnroHyr o SUB TABLE 
. SET UP TABLES 


[xqt aus 

TABLE(NI=3,NJ=1 , ITYPE=0) : SUBS PARA 1 1 
J=l: 2 4 1 . nsubs nnpe iopt 

« 

. NO. OF ELEMENTS IN SUBSTR. 

TABLE(NI=2,NJ=1 , ITYPE=0) : SUBS ELEM 1 1 

J=l: 4 4 .no. of elements per sub. (for n subs) 

. SUBSTRUCTURE INTERFACES 

. syntax: 

. J=m: jl si j2 s2 sn=l,no. of interface nodes 

TABLE(NI=4,NJ=3.ITYPE=0): SUBS SUBI 1 1 
M: 3 112 
J=2: 6 12 2 
J=3: 9 13 2 

*end 

*call SUBJABLE 
*add preserye 


. DEFINE SUBSTRUCTURE -1 

[xqt tab 
start 9, 4, 5,6 
mate:! 10000 0.3 
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ORIGINAL’ PAGE IS 
oniine=0 OF POOR QUAUTX 

jloc: format=l 

1 0.0 0.0 0.0 15.0 0.0 0.0 3 i 3 
3 0.0 11.5 0.0 15.0 11.5 0.0 
sa(l):i 0.1 
con=i 

zero 3:1,9 
[xqt eld 
onl ine=0 
e41 

1 2 5 4 1 2 2 
online 2 ! 

[xqt topo 
[xqt e 
[xqt eks 
stop 

★call end_sub 

a 

. DEFINE SUBSTRUCTURE -2 


[xqt tab 
start 9, 4, 5,6 
mate: 1 10000 0.3 
jloc: format=l 
online=0 

1 15. 0. 0. 15. 11.5 0. 3 1 

4 22.5 0.0 0.0 30.0 0.0 0.0 2 1 3 

2 22.5 11.5 0.0 30.0 11.5 0.0 
sad) : 1 0.1 

000=1 

zero 3:1,9 
zero 1:5, 9, 2 
zero 2:9 
[xqt eld 
online=0 
e41: 14 6 2 
2 6 8 3 
4 5 7 6 1 1 2 
online=l 
[xqt topo 
[xqt e 
[xqt eks 
stop 

*caii end_sub 
★call giob_dat 


[xqt k 
[xqt inv 
[xqt aus 
alpha 
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case titles 
r external load 
sysvec 

applied forces 
case 1 
1=1 

J«l:2.875 
J=4:5.75 
J=7:2.875 
[xqt ssoi 
[xqt gsf 
[xqt dcu 
[xqt osf 
[xqt vprt 
tprint stat disp 
tprint stat reac 
[xqt exit 
endinput 


ORIGINAL PAGE IS. 
OF POOR ouAura 



ORIGINAL PAGE IS 
OF POOR QUALITY 

OUTPUT FOR EIGHT-ELEMENT RECTANGULAR PANEL 


<Cl> PUTjrsess sqe , Commn t > 

<CL) $root f LC001 ,CuGGQi>*set echo=off 

urEN, Ldi: 1, Files SSi.LOl , Attr: nsw, Block I/O 
** BEGIN TAB ** DATA SPACE* 200000 WORDS 

idUU JGINiS. 

0 ACTIVE JOINT MOTION CGHPGNENTS* 1 2 3 4 5 5 
Oils READ ERROR, NU , L , XCRE , I ERR , 1 0P= 1 Q 196333 

NAME* MASK STAB 2 5 

0*** ERRORS IN INPUT PREVENT CALCULATION OF QJ(3,3,JT) 

EXIT T AB CPUTIME* 7.0 I/0(DIR,SUF)= 0 0 

** BEGIN AUS ** DATA SPACE* 200000 WORDS 

TABL COMPLETED. 

TABL COMPLETED. 

TABL COMPLETED. 

EXIT AUS CPUTIME* 5.3 i/0(DIR,BUF)= 0 0 

** BEGIN AUS irk DATA SPACE* 200000 WORDS 

MACR COMPLETED. 

MACR COMPLETED. 

MACR COMPLETED. 

EXIT AUS CPUTIME* 1.3 I/0(DIR,3UF5* 0 0 

3 JOINTS. 

UHLiivc, JuiNT MOT i UN COMPONENTS* 12 3 
EXIT TAB CPUTIME* 47.6 I/0(DIR,SUF)= 0 0 

** BEGIN ELD ** DATA SPACE* 200000 WORDS 

EXIT ELD CPUTIME* 7,4 I/0(DIR,BUr)= 0 0 

** BEGIN TOPS ** DATA SPACE* 200000 WORDS 

■sDM! OPEN, Lei: 25, File: NS.l 25 , Attr: scratch, Slock I/O 
<DM) OPEN, Ldi: 25, File: N3.L2S , Attr: scratch, Slock I/O 
NO, OF 4-NODE ELEMENTS* 4 
0TGTAL NO. OF ELEMENTS* 4 
OrAXCQN, MAXSUB, ILMAX* 1872 1400 52 


0KSIZE,NR5,LR5= 

13 1 

336 

OMAXCON, MAXSUB, ILM 

AX* 1864 1400 

52 

3SI2E INDEX* 15, 

I Cl, ICG* 35 

33, NR4* 1 

<DM) CLOSE, Ldi: 25 

, File: NS.L25 


<DM) CLOSE, Ldi: 26 

, i- lie: NS.L26 


EXIT TOPO CPUTIME* 

7.3 I/O (DI R , SUF) 

= 0 0 

"rk BEGIN E k* 

DATA SPAuE* 20000 

iO WORDS 

"= 0.100GGE-13-0. 

10000E-02 Q.10000E- 

■04 u . IflOyOE-O* 

0.20000Et02 0. 

iuuGQE-jj d.luuuOt- 

■03 0.10000E-03 

ERROR LEVELS* 2 2 0 

2 2 2 2 2 



L, VOL OR 
3UP AREA SLM WEIGHT 
1 C.’72500Er03 0.0 000 ODE’- 00 


STRUCTURAL NCN— STRUCTURAL 

WEIGHT 

ii i y v w L : v iTL w 


0 .OOOOOOErOO O.OOSOOOErOO 
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i 


ORIGINAL PAGE 13 
OF POOR QUALITY 


i PfL 




./Ui UiK,iUr }- 


EKS ** DATA SPACE- 20G0U0 WORDS 
1 , 3 I/O (DIR, BUF; - G 




LAii uhi LfUi* 

Uj\H ? i 

** BEGIN AU3 ** DATA SPACE 5 2G0090 WORDS 

•••ACS COMPLETED. 
MACR COMPLETED. 
MACR COMPLETED# 
MACS COMPLETED, 

OKAY 1.2 
QKAYY 1.4 
njsSCt = 43 
nssctli] = 1 


MhCR CQHP‘.£TED. 


EXIT AUS CPUTIME = 3.3 I/0(DIR,BUF> 


;?•, • ; .r. . . ',i - q 

V w W / 1 r I V V I W I C 

OKAY 2 l 2 
QKAYY 3 3 3 




13.300000 


2S.3 I/0(Di?.,3UF)=' 


TATi 


.3 I/0( 

CLCilN i Uf U ** L'H ; « orH^L- iUUuUU WuiWi 

•:c«> GPE?4 f Ldi : 25, File: NS.L25 , At tr ; scratch, Block 
< DM ) OPEN, Ld 1 : 26, File! NS . u2b 
NO. OF 4-NODE ELEfiENTS= 4 
CTjTAl NO. OF ELEMENT 4 
urkX C uN , MAX cub, * ltlAX 5 1 8 .» L 


4U’j 


3KS:ZE,NR5,lR5= 


0MAXCQN, MAXSU8 , I LHAX= 1364 14 

0SIZE INDEX 5 15, 101, IC2= 

<3M) CLOSE, cdi: 25, Fils: NS.L25 
vDn) CLOSE, Lei: 26, File: NS.l26 


83b 

32, NR4= 


EXIT TOPO CPUTIME 5 10.8 I/Q(DIR,EliF} = 0 

•At BEGIN E ** DATA SPACE 5 23G003 WORDS 
T= 0 .10G00E-13-0 .10C00E-02 G.10COOE-04 O.iGOOGE-04 
0.2 jD00EtG2 0.10000E-03 O.1QO80E-Q3 0.10GOGE-03 
ERROR LEVELS 5 22 0 22222 


L f VOL OR 
Tr?E GROUP AREA Sui 
E41 1 0.172500 EtG 3 O.OOOGOOEtOO 0 


STRUCTURAL NCM-STRUCTURAl 
wEIu.HT MEiliHi 

*J u v v u Lr b U 


J ; riL 


0 , QQCGC OEAOO O.OOOOOGEtOO ■ 
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TOTAL 4-NODE? 


CPUTIME 3 


** BEGIN EKo kk DATA SPACE 3 


250CGE+03 

10.7 i/0(DIR,3UF)= 0 

200000 WORDS 


COMPLETED 
EK3 CPUTIME 3 


1/ 0 { Dl K i dJF; - 


** BEGIN AU3 ** 


DATA SPACE 3 200000 'WORDS 
MACR COMPLETED . 
MACR COMPLETED. 
MACR COMPLETED. 
•MACR COMPLETED . 


OKAY 1.2 
GKAYY 1.4 

j zzrt - *2 

•’SfctiO] 3 1 

EXIT AUS CPUTIME 3 
<CL> lvsO t : G1 5 . 8 
** BEGIN AUS rk 


MACR COMPLETED. 
12.1 I/O (DIR, 5UF) = 0 

SECTLE 1 ] 10,000000 

DATA SPACE 3 200000 wORDS 





MACR 

COMPLETED 




MACR COMPLETED 




MVACR COMPLETED 




MACR COMPLETED 

(CD 

iv:0 

t : I 

GJNT[1] 

j. 

<CL> 

lv:0 

t : I 

GJNT[2] 

2 

<CL> 

ly : U 

t : I 

GJNTC33 

3 

(CL) 

iv:G 

tsl 

GJNT[ 4] 

4 

<CL> 

lv:0 

t i ; 

GJNTi.5] 

j 

/ r- \ 

iv :0 

til 

GJNT[t] 

6 

(CD 

1 y : G 

tsl 

GJNT[7] 

/ 

(CD 

lv : 0 

til 

3JT4T[ 6] 

8 

/■ 

lv :U 

t : I 

GJNTE3] 

3 

V-L / 

D;0 

til 

SJNT l 10 j 


pi V 

) , ;* 5 

* t T 

^ . 4 4 , 


w w i 


i ; i 

tuN. L-i J 


\t,L> 

1m:0 


QJNT[12] 

3 

\CL) 

iv :0 

t : 1 

GJNT[13] 

i-J 

(CD 

lv sO 

tsl 

GJNT[14] 

11 

‘-.vL/ 

lv : 0 

t i i 

GJNT[15] 

* ^ 
mL 

\ L w / 

lv sO 

ti! 

GJNT: le‘ 

13 

<CL) 

lv :u 

t:I 

SJNTil?] 

14 

\ L-^/ 

iv ;u 

t i i 

GJNTEIS] 

i - 

EXIT 

AUS 

f-p *:TTyr- 

wf A ; i. 

130.1 1/0 (DIR, 

3UF)= 


** BEGIN AUS ** 


UfiiH opHlt 3 iUlUUO wukUc 

DEF5 COMPLETED. 
DEFI COMPLETED. 
TA3L CuiFLETED . 
MACR COMPLETED. 
TA8L COMPLETED. 
TA3L COMPLETED. 
TA8L CCMP-ETED. 
TABL DEPLETED . 
TA6L COMPLETED. 

Lt'l 


DMPi i 


ORIGINAL PAGE IS 
3 OF POOR QUALITY 
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EXIT AUS CPUTIME 8 53,3 I/Q(DIR,SUF)= Q 0 

** BEGIN AUS ** DATA SPACE 8 200000 WORDS' 

DEFI COMPLETED. 

DEFI COMPLETED. 

TABL COMPLETED. 

EXIT AUS CPUTIME 8 13. 6 I/OCDIR.BUF) 8 0 0 

START 13 4,5,6 

EXIT TAB CF'uTIME 8 11.5 I/O f DI R ,SUF) = 0 0 

** BEGIN AUS ** DATA SPACE 8 200000 WORDS 

MACR COMPLETED. 

$ of elements in substructure 1-4 
NJ of DEF dataset is 56 

MACR CCMPLETED. 

* of elements of substructure 2-4 

MACR COMPLETED. 

TABL COMPLETED. 

MACR COMPLETED. 

TABL COMPLETED . 

MACR COMPLETED. 

TABL COMPLETED. 

MACR COMPLETED. 

TABL COMPLETED. 

EXIT AUS CPUTIME 8 63.3 I/0(DIR,BUF5 8 0 0 

** BEGIN ELD ** DATA SPACE 8 200000 WORDS 

EXIT ELD CPUTIME 8 16. S I/O ( 01 R , 8UF) 8 0 0 

** BEGIN AUS irk. DATA SPACE 8 200000 WORDS 

TABL COMPLETED. 

EXIT AUS CPUTIME 8 4.0 I/0(DIR,BUF)= 0 0 

** BEGIN TOFO kx DATA SPACE 8 200000 WORDS 
<DM> OPEN, Ldi: 25, File: NS.L25 , Attr: scratch. Sleek 1/C 
<DM> OPEN, Ldi: 26, File: NS.L26 , Attr: scratch, Block I/O 
NO. OP 4-NODE ELEMENTS 8 3 
OTOTAL NO. OF ELEMENTS 8 S 


Hi V r- ! M w**| in 

U! t-IAC.Ul'l , f!HA3U0. 

ILMAX 8 

1372 

1400 

52 

QKSIZE,NR5,LR5 8 

10 


1 

836 

OMAXCON, MAXSU3, 

I lMAX 8 

1364 

1400 

cr 

wCl 

OSIZE INDEX 8 

23, I Cl 

, IC2= 

221 

70, NR4= 


<DM> CLOSE, Ldi: 25, File: NS.L25 

tOM) CLOSE, Ldi: 26, File: N3.L26 

EXIT TOPS CPUTIME 8 18.0 i/0(DIR,BUF)= 0 0 

irk iEGIN AUS krk DATA SPACE 8 200000 WORDS 

DEFI COMPLETED. 

DEFI COMPLETED. 

UNIO CCMPLETED. 

EXIT AUS CPUTIME 8 7.3 I/0(DIR,8UF5= 0 0 

** BEGIN K ** DATA SPACE 8 2COOQO WORDS 

EXIT K CPUTIME 8 10.3 I/0(DIR,8UF)= 0 0 

*k BEGIN INV ** DATA SPACE 8 200000 WORDS 

GNSING.NNEG 8 0 0 

EXIT INV CPUTIME 8 6.3 I/S(DIR,3Ur) 8 0 0 

** BEGIN AUS ** DATA SPACE 8 200000 WORDS 


ORIGIN All PAGE IS 
OE POOR QUALITY 


3YSV CCMPLETED. 
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EXIT AUS CPUTIME- 4.8 I/0(DIR,3UF)= 
bEliiN S50L ** DATA SPACE- iiOuOOO WORD* 


<CH> OPEN, Ldi: 

26, File: 

N3.L2S , Attr: 

scratch, 

0 CASE 

F*U 


U*KU 

ERR 

1 

Q.345GQGQE+G0 0, 

,3450001E+G0 -0, 

, 863S353E 

(DM) Ci 

-GSE, Ldi: 

26, File; 

: NS.L26 



EXIT 3SGL C?U7IME= 16,3 I/O ( DIR, BUF) = 0 0 

** BEGIN DCli ** DATA SPACE= 200000 WORDS 

0 

EXIT DCU CPUTIME 5 0.2 I/0(DIR,8UF)= 0 0 

** BEGIN VPRT ** DATA SPACE= 200QG0 WORDS 
1 STATIC DISPLACEMENTS, 
external loao 


0 JOINT 1 

i, 

3 

i 0.30GE-01 

"U .o*4-jL~uc: 

v ivuwL Uu* 

*i MMcr rsi 
i •; .Uwt‘Ui 

rt hjcr rt-n 
“U . OtJE.~Vfc 

G.iiGGE+uG* 

3 Q.15QE-Q1 

-0.345E-02 

G.QGQE+00* 

4 0.300E-01 

"0 .173E-02 

O.OOOE+OO* 

5 9.225E-01 

-0.173E-02 

0 .OOOE+OO* 

6 0.150E-01 

-0.173E-02 

O.OOOE+OO* 

7 0.300E-01 

-0.753E-07 

O.OOOE+OO* 

3 0.225E-01 

-8.521E-0? 

O.OOOE+OO* 

3 0.150E-B1 

-0.236E-07 

O.OGQE+03* 

10 0.75GE-G2 

-G.345E-02 

O.OOOE+OO* 

11 O.OOOE+OO* 

-0.345E-02 

O.OOOE+OO* 

12 0.750E-02 

-0.173E-02 

O.OOOE+OO* 

13 G.GDGEtOQ* 

-Q.173E-02 

0.000E+00* 

14 0.750E-02 

-0.120E-07 

O.OOOE+OO* 

15 O.QOOE+OU* 

0 .OODE+Oo* 

0.000E+00* 

1 STATIC REACTIONS, 

FORCE ERRORS. 

external load 

.* rr-.r vrr * 



w'UUiSt: 1 

£ 

0 

1 — 0 .113E-05 

-G.283E-QS 

Q.jQuE+00* 

2 Q.62SE-06 

0.113E-06 

Q.000E+GG* 

3 -C.121E-05 

0.243E-06 

0.000E+0C* 

4 3.143E-05 

0.586E-07 

G.QOGE+00* 

5 Q . 40GE-GO 

-0.143E-06 

O.OOOE+OO* 

A J-w * A f— 

b -U.b^iL-U/ 

-0.103E-07 

0.000E+00* 

r. a r 

i ~U ,4;'/'t”Ub 

0.533E-06 

O.OOOE+OO* 

3 -0.212E-3S 

-0.334E-36 

O.OOOE+OO* 

3 -0.9S3E-Q7 

-G.S30E-07 

O.OOOE+OO* 

10 0.537E-0S 

-0 .335E-06 

0 i vuOw+UO* 

ii -u.283Et01* 

-0.634E-07 

o .U cOEtOO* 

12 0.563E-SS 

A ~ - 

U.£?UL-lib 

O.OOOE+OO* 

13 -Q.STSE+Qi* 

0 ,lbr E“06 

O.GOOE+uO* 

14 -0.233E-07 

U ,67it~07 

Q, OOOE+OO* 

15 -0.287E+01* 

0.178E-05* 

O.OOOEtOO* 

EXIT vPRT CPUTIME 1 

= 4 . L I/u( 

DIR,3UF)= 

<DM> CLOSE, wGi : 

1, File: SSI 
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!!! link with /usr/ul/elp/lifa/nslib.a !!!!! 

Force LliD of NP ident ME 

Test program for parallel PLU decomposition using Force versions 
of *Linpack‘ routines. 

This program calls a (single stream) routine to read in a matrix 
and vector, then SGEFA is called to factor the matrix into a PLU form 
then SGESL is called to solve the system. The computation is timed, 
and then a (single stream) routine is called to output the results. 

Shared INTEGER N, IPVT(500), INFO 
Shared DOUBLE PRECISION A( 500, 500) ,B(500) 
shared INTEGER ITM1, ITM2, ITM3 
shared REAL DELTA1 , DELTA2, RATE1, RATE2 
Private integer amax, numjnt,dof,ierr, rmax, bw 
Private integer rowptr(lOOO), rowlen(lOOQ) 

Private double precision amat(5QQOO),x(lQ0O) 

Private integer i, j 
End declarations 

Barrier 
amax = 10000 
rmax = 1000 

call getrow( ama t , amax , rowp tr , r owlen , rmax , 
'/usr/ul/venkatesh/panel/MAS.LQl' , 
x,b,numjnt,dof,bw,ierr) 
if (ierr.ne.O) write!*, 5) ierr 
format('error return code = ',i2) 

N = numjnt*dof 

print*, 'starts to print matrix' 
do 1 i * 1,N 

write(*,2) i,rowptr(i),rowlen(i),amat(rowptr(i)) 
continue 

format('row ',i3,' start ',i6,'lth= ',i3,'diag= ',gl5.5) 

...Code to transfer from vector to symetric array, a 
print*, 'starts symmetric array' 
do 20 irow = 1,N 
ip = rowptr(irow) 
ilen = rowlen(irow) 
ilen = (irow-1) + rowlen(irow) 
do 10 icol = irow, ilen 
a(irow,icol) = amat(ip) 
ip = ip + 1 
continue 

zero out the remainder of the array 

do 15 icol = ilen+1, N 
a(irow,icol) = 0. 
continue 
continue 

print*, y o k a y 3' 
do 40 icol = 1, N 
do 30 irow = icol, N 
a( irow, icol) = a(icol,irow) 
continue 
continue 


C CALL PRNTB2CA, B, N) 

CALL CFrtic(ITMl) 

End barrier 

Forcecall 3GEFA(A, 500, N, I PVT, INFO) 

IF (INFO .NE. 0) THEN 
C matrix is singular... 

if (me .eq. 1) print *, 'matrix is singular' 

Join 
END IF 

print*, 'o k a y 4' 

IF (ME .EQ. 1) CALL CFrtic(I7M2) 

Forcecall SGESL(A, 500, N, IPVT, B, 0) 

Barrier 

print*, 'o k a y 5' 

CALL CFrtic(ITM3) 

CALL PRINTB(A, B, N) 

DELTA1 = (I7W2 - ITM1)/50.0 
DELTA2 = (ITM3 - ITM2)/50.0 
OPEN( 8 , F I LE=*/u sr/ul/v enkatesh/time.dat',CPU=l) 

IF (DELTA1 .NE. 0) THEN 
RATE1 = ( 2 . Q*N*N*N/3 . 0 ) *1 . OE-fi/DELTAl 
WRITE (8, 400) NP, N, ME, DELTA1, RATE1 
ELSE 

WRITE (8, 400) NP, N, ME, DELTA1 , 0.0 
WRITE (8, 402) 

ENDIF 

IF (DELTA2 .NE. 0) THEN 
RATE2 = (N*N)*1 .0E-6/DELTA2 
WRITE (8, 401) ME, DELTA2, RATE2 
ELSE 

WRITE (3, 401) ME, DELTA2, 0.0 
WRITE (8, 402) 

ENDIF 

print*, 'o k a y S' 

400 FORMATS Number of processes =', 13,' Matrix order =', 14,/ 

+ ' matrix decomposition: ',F12.4,' seconds (',F5.3,' MFLOPS)') 

401 FCRMAT( 

t ' Fwd/back subst: ',F12.4,' seconds (',F5.3,' MFLOPS)') 

402 FORMAT!' No elapsed time measured.') 

End barrier 

Join 

END 
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c 

c 
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c 
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20 
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40 


Forcesub SGESL(A, LDA, N, I PVT, 8, JOB) of NP ident me 
3GESL; Force parallel version of UNPACK' s SGESL which solves 
Ax=b, using (column orientated) foward/backward subst, 

A call to SGEFA must precede the call to SGESL. 

SGEFA provides the PLU factorization of A. 

A: the working matrix 

LDA: the declared dimension of A 

N: the working dimension of A 

IPVT: the pivot array 

B: the input/output vector (solved in place) 

B is the only argument altered by SGESL. 

JOB: job selector: job=0: solves Ax=b 

job=l: solves Transpose(A)x=b, not implemented 
Ail argument MUST be declared ’Shared* in the calling module! 
integer LDA, N, IPVT(LCA), JOB 
double precision A(LDA,LDA), B(LDA) 

Private integer I, K, L 
Shared real TEN 
End declarations 

first solve: L*Y = B 
do 20 K = 1, N-l 
Barrier 

L '= IPVT(K) 

TEM = B(L)' 

B(L) = B(K) 

8(K) = TEM 
End banier 

Presched- do 22 I = K+l, N 
B(I) = B(I) + TEM * A(I,K) 

End presched do 
continue 

then solve: U*X = Y 
do 40 K = N, 1, -1 
Barrier 

8(K) = B(K)/A(K,K) 

TEM = -B(K) 

End barrier 

Presched do 42 1=1, K-l 
B(I) = B(I) + TEM * A(I,K) 

End presched do 
continue 
Barrier 
End barrier _ 

RETURN 

END 


Forcesub SGEFA(A, LDA, N, I PVT, INFO) of NP ident HE 
INTEGER LDA, N, IPVT(LOA), INFO 
DOUBLE PRECISION A(LDA, LDA) 

C 

Private REAL ABSMAX 
Private REAL T 

Private INTEGER KK, I, MAXI, J, KP1 
Shared INTEGER IALL 
Shared REAL PIV, ALLMAX 
Shared logical ILOCK 
End declarations 

Barrier 
INFO = 0 
ALLMAX = 0.0 
End barrier 

C.......For each row of the matrix 

DO 1000 KK = 1, N-l 
IF( INFO.NE.O) GOTO 2000 
C....... Reset local maxima 

ABSMAX = 0.3 

C.......Find the pivot row 

> 

Presched DO 100 I = KK, N 
T = ABS( A(I, KK) ) 

IF (ABSMAX .LT. T) THEN 
MXI = I 
ABSMAX = T 
ENDIF 

100 End Presched DO 

C Update the Shared Information if necessary 

Critical ILOCK 

IF (ALLMAX .LT. ABSMAX) THEN 
IALL = MAXI 
ALLMAX = ABSMAX 
ENDIF 

End critical 
Barrier 

IF (ALLMAX .EQ. 0.0) THEN 
INFO = KK 
IPVT(KK)=KK 
ELSE 

iPVT(KK)=IALL 

ENDIF 

End Barrier 
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If the matrix is singular then pass the information 


IF (INFO .EQ. KK) GOTO 1000 
MAXI = I ALL 

IF (MAXI .NE. KK) THEN 

C Swap rows if necessary 

Presched DO 110 J = KK, N 
TEMP = A(MAXI , J) 

A(MAXI, J) = A(KK, J) 
A(KK, J) = TEMP 


110 End Presched DO 

ENDIF 

C.......Self schedule row reductions 


KP1 = KK + 1 
Barrier 

PIV = -1.0/A(KK, KK) 

End Barrier 

Selfsched Do 130 I s KP1, N 
TEMP = PIV*A(I, KK) 

A(!,KK) = TEMP 
DO 120 J = KP1, N 

A(I, J) = A(I, J) + TEMP*A(KK, J) 

120 CONTINUE 

130 End Selfsched DO 

Barrier 
ALLMAX=0 
End barrier 

1000 CONTINUE 

2000 RETURN 

END 

**»** *AAAAAAAAAAAAAAAA A^ * *** A*AAAAAAAAAAAAAAAAAAAAAAAAA *THHHt- A - A* - A * AA * AAAAAAAAA 
AA AAAAAAAAAAAAAAA AAAArAAAAA- A i/o routines (replacable) AAAAAAAAAA ' AAAA A AAA A AAAAAAA 
AAAAAAAAAAAA A**** A ** A* A*A** AAAAAAAAAA*AA **A********rt********** AAA*AAA *** * **- A A* 


subroutine printb(a, b, n) 

C this routine prints 8(N) & A(N,N) 

C..... .where B(N) is the solution vector, and A(N,N) is the working matrix. 

C the user may replace this single stream subroutine. 

double precision b(i), a(500,500) 
open(8,FILE=*/usr/ul/venkatesh/panei/venk .dat 1 ,CPU=1) 
rewind (8) 
do 10 i=l,n 

10 write(3,90) ^b( ' , i , ' )=' ,b( i ) 

close (3) 

90 format(a,i3,a,el4.7) 

return 
end 



subroutine prntb2(a, b, n) 
this routine prints 3(N) & A(N,N) 

where 3(N) is the solution vector, and A(N,N) is the working matrix, 
the user may replace this single stream subroutine, 
double precision b(l) , a(5G0,500) 
open(8,FIL£= , /usr/ui/venkatesiv'panei/venk.pred3t',CPlj=15 
rewind (8) 
do 10 i=l,n 

write(8,30) 'b(',i/)=',b(i),(a(i,j), j=l,n) 
ciose (8) 

format(a,i3,a,300el4.7) 

return 

end 


OUTPUT FOR RECTANGULAR PANEL WITH HOLE 


(Cl> PUT_message,Coiwnnt> 

<CL> $root,L0001,C00001>*set echo=off 

<DM> OPEN, Ldi: 1, File: PW1.L01 , Attr: new, Block I/O 

** BEGIN TAB ** DATA SPACE= 200000 WORDS 

1000 JOINTS. 

0ACTIVE JOINT MOTION COMPONENTS 3 1 2 3 4 5 6 

0LIB READ ERROR, NU,l,KORE,iERR,iOP= 1 0 158988 -1 11 

NtfiE 3 hVASK BTAB 2 5 

0*** ERRORS IN INPUT PREVENT CALCULATION OF QJ(3,3,JT) 

EXIT TAB CPUTIME 3 7.5 I/0(D!R,BUF) 3 0 0 

** BEGIN AUS ** DATA SPACE 3 200000 WORDS 

TABL COMPLETED. 

TABL COMPLETED. 

TABL COMPLETED. 

EXIT AUS CPUTIME 3 5.7 I/0(DIR,BUF) 3 0 0 

** BEGIN AUS ** DATA SPACE 3 200000 WORDS 

MCR COMPLETED. 

WCR COMPLETED. 

MACR COMPLETED. 

EXIT AUS CPUTIME 3 1.8 I/0(DIR,BUF) 3 0 0 

24 JOINTS. 

OACTIVE JOINT MOTION COMPONENTS 3 1 2 3 
EXIT TAB CPUTIME 3 50.0 I/0(DIR,BUF)= 0 0 

LREC 3 1280 

** BEGIN ELD ** DATA SPACE 3 200000 WORDS 

EXIT ELD CPUTIME 3 8.8 I/0(DIR,BUF) 3 0 0 

** BEGIN TOPO ** DATA SPACE 3 200000 WORDS 

<DM> OPEN, Ldi: 25, File: NS.L25 , Attr: scratch, Block I/O 

<DM> OPEN, Ldi: 26, File: NS.L26 , Attr: scratch, Block I/O 

NO. OF 4-NODE ELEMENTS 3 16 


OTOTAL NO. OF ELEMENTS 3 

16 



OMAXCON, MAXSUB, ILMftX 3 

1872 

1400 

52 

0KSIZE,NR5,LR5= 28 


1 

896 

OMAXCON, MAXSUB, ILMAX 3 

1864 

1400 

52 

OSIZE INDEX 3 66, IC1 

, IC2= 

967 

i 3 

II 


<DM> CLOSE, Ldi: 25, File: NS.L25 
<DM> CLOSE, Ldi: 26, File: NS.L26 
EXIT TOPO CPUTIME 3 7.2 I/0(DIR,BUF)= 0 0 

** BEGIN E ** DATA SPACE 3 200000 WORDS 
T= 0 .10000E-19-0 .10000E-02 0.10000E-04 0.10000E-04 
0.20000E+02 0.10000E-03 0.10000E-03 0.10GQOE-03 
ERROR LEVELS 3 2 2 0 2 2 2 2 2 

0 

0 

L, VOL OR STRUCTURAL NON-STRUCTURAL 

TYPE GROUP AREA SUM WEIGHT WEIGHT 

E41 1 0.131 71 6E+02 O.OOQOOOE+OO D.OOOOOOE+OO 



O.OOOOOOE+OO O.OOOOOOE+OO 


0 

TOTAL 


TOTAL 4-NODE? 0.1317158E+Q2 

EXIT E CPUTIME= 6.7 I/0(DIR,BUF)= 0 0 

** BEGIN EKS ** DATA SPACE= 200000 WORDS 

E41 COMPLETED 

EXIT EKS CPUTIME* 6.7 I/0(DIR,BUF)= 0 0 

OKAY 1 

** BEGIN AUS ** DATA SPACE= 200000 WORDS 

MACR COMPLETED. 

MCR COMPLETED. 

NO COMPLETED. 

MACR COMPLETED. 

OKAY 1.2 
OKAYY 1.4 
njsect = 43 
nsect[l] = 1 

MACR COMPLETED. 

EXIT AUS CPUTIME* 9.0 I/0(DIR,BUF)= 0 0 

<CL> lv:0 t:G15.8 SECTLfl] 10.000000 
OKAY 2 2 2 
OKAYY 3 3 3 

EXIT TAB CPUTIME* 42.4 I/0(DIR,BUF)* 0 0 

** BEGIN ELD ** DATA SPACE= 200000 WORDS 

EXIT ELD CPUTIME* 13.8 I/0(DIR,BUF)= 0 0 

** BEGIN TOPO ** DATA SPACE* 200000 WORDS 

<DM> OPEN, Ldi: 25, File: NS.L25 , Attr: scratch, Block I/O 

<DM> OPEN, Ldi: 26, File: NS.L26 , Attr: scratch, Block I/O 

NO. OF 4-NODE ELEMENTS* 44 


OTQTAL NO. OF ELEMENTS* 

44 



OMAXCON, MAXSUB, IL+AX* 

1872 

1400 

52 

OKSIZE ,NR5 , LR5= 52 


2 

896 

OMAXCON, MAXSUB, ILMAX* 

1863 

1400 

52 

OSIZE INDEX* 210, IC1, 

, IC2* 

4905 

702, 1 


<DM> CLOSE, Ldi: 25, File: NS.L25 
<DM> CLOSE, Ldi: 26, File: NS.L26 • 

EXIT TOPO CPUTIME* 14.9 I/0(DIR,BUF)= 0 0 

** BEGIN E ** DATA SPACE* 200000 WORDS 
T= 0 .10000E-19-0 .10000E-02 0.10000E-04 0.10000E-04 
0.20000E+02 0.10000E-03 0.10000E-03 0.10000E-03 
ERROR LEVELS* 2 2 0 2 2 2 2 2 

0 

0 

L, VOL OR STRUCTURAL NON-STRUCTURAL 

TYPE GROUP AREA SUM WEIGHT HEIGHT 

E41 1 0.329000E+03 O.OOOOOOE+OO O.OOOOOOE+OO 

0 

TOTAL 


O.OOOOOOE+OO O.OOOOOOE+OO 



TOTAL 4-NODE# 0.3290000E+03 

EXIT E CPUTIME 2 11.7 I/0(DIR,BUF)= 0 0 

** BEGIN EKS ** DATA SPACE 2 200000 WORDS 

E41 COMPLETED 

EXIT EKS CPlfTIME= 15.9 I/0(DIR,BUF) 2 0 0 

OKAY 1 

** BEGIN AUS ** DATA SPACE 2 200000 WORDS 

mCR COMPLETED. 
fttCR COMPLETED. 

MACR COMPLETED. 

MACR COMPLETED. 

OKAY 1.2 
OKAYY 1.4 
njsect = 43 
nsect(2] = 1 

WCR COMPLETED. 

EXIT AUS CPUTIME 2 12.1 I/0(DIR,BUF)= 0 0 

<CL> 1m: 0 t:G15.8 SECTLtlJ 10.000000 
OKAY 2 2 2 
OKAYY 3 3 3 

** BEGIN AUS ** DATA SPACE 2 200000 WORDS 

MCR COMPLETED. 

MACR COMPLETED. 
fttCR COMPLETED. 

M4CR COMPLETED. 

EXIT AUS CPUTIME 2 954.2 I/0(DIR,BUF)= 0 0 

** BEGIN AUS ** DATA SPACE 2 200000 WORDS 

TABL COMPLETED. 

TABL COMPLETED. 

TABL COMPLETED. 

EXIT AUS CPUTIME 2 214.2 I/0(DIR,BUF)= 0 0 

** BEGIN AUS ** DATA SPACE 2 200000 WORDS 

DEFI COMPLETED. 

DEFI COMPLETED. 

TABL COMPLETED. 

EXIT AUS CPUTIME 2 25.3 I/0(DIR,BUF) 2 0 0 

START 78 4,5,8 

78 JOINTS. 

OACTIVE JOINT MOTION COMPONENTS 2 1 2 3 
EXIT TAB CPUTIME 2 10.7 I/0(DIR,BUF) 2 0 0 

** 8EGIN AUS ** DATA SPACE 2 200000 WORDS 

MCR COMPLETED. 

* of elements in substructure 1-16 
NJ of DEF dataset is 30 

MACR COMPLETED. 

# of elements of substructure 2-44 

** BEGIN ELD ** DATA SPACE 2 200000 WORDS 

EXIT ELD CPUTIME 2 20.7 I/0{DIR,8UF)= 0 0 

** BEGIN AUS ** DATA SPACE 2 200000 WORDS 

TABL COMPLETED. 
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0 


EXIT AUS CPUTIME= 4.2 I/O (DI R , BUF) = 0 

** BEGIN TOPO ** DATA SPACE 1 2 200000 WORDS 
<DM> OPEN, Ldi: 25, File: NS.L25 , Attr: scratch, 81ock I/O 
<DM> OPEN, Ldi: 26, File: NS.L26 , Attr: scratch, Block I/O 
NO. OF 4-NODE ELEMENTS 2 60 
OTQTAL NO. OF ELEMENTS 2 60 


OMAXCCN, MAXSUB , 

, ILMAX 2 1872 

1400 

52 

0KSI2E,NR5,LR5= 

52 

3 

896 

OMAXCCN, MAXSUB, 

, ILMAX 2 1863 

1400 

52 

0SI2E INDEX 2 

210, IC1, IC2= 

5923 

872, NR4= 4 


<DM> CLOSE, Ldi: 25, File: NS.L25 

<DM> CLOSE, Ldi: 26, File: NS.L26 

EXIT TOPO CPUTIME 2 21.2 I/0(DIR,BUF)= 0 0 

** BEGIN AUS ** DATA SPACE 2 200000 WORDS 

DEFI COMPLETED. 

DEFI COMPLETED. 

OHIO COMPLETED. 

EXIT AUS CPUTIME 2 15.1 I/0(DIR,BUF)= 0 0 

** BEGIN K ** DATA SPACE 2 200000 WORDS 

EXIT K CPUTIME 2 28.0 I/0(DIR,BUF) 2 0 0 

** BEGIN INV ** DATA SPACE 2 200000 WORDS 

QNSING,NNEG 2 0 0 

EXIT IW CPUTIME 2 20.7 I/0(DIR,BUF) 2 0 0 

** BEGIN AUS ** DATA SPACE 2 200000 WORDS 

SYSV COMPLETED. 

EXIT AUS CPUTIME 2 3.5 I/0(DIR,BUF)= 0 0 

** BEGIN SSOL ** DATA SPACE 2 200000 WORDS 

<DM> OPEN, Ldi: 26, File: NS.L26 , Attr: scratch, Block I/O 

0 CASE F*U U*KU ERR 


1 0.3532554E+00 0.3532554E+00 -Q.1687296E-06 

<DM> CLOSE, Ldi: 26, File: NS.L26 
EXIT SSOL CPUTIME 2 18.1 I/0(DIR,BUF)= 0 0 

** BEGIN DCU ** DATA SPACE 2 200000 WORDS 

0 


EXIT DCU CPUTIME 2 0.2 I/0(DIR,BUF)= 0 
** BEGIN VPRT ** DATA SPACE 2 200000 WORDS 


OLIB READ ERROR, NU,L,KORE,IERR,IOP= 
WME= CASE TITL 11263747405 

1 STATIC DISPLACEMENTS. 


1 


QJGINT 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


1 

0.123E-01 

0.134E-01 

0.152E-Q1 

0.174E-01 

0.134E-01 

0.174E-Q1 

0.152E-01 

0.134E-01 

0.123E-01 

0.136E-01 

0.154E-01 


-0.173E-02 

-0.107E-02 

-0.716E-03 

-0.115E-02 

-0.173E-02 

-0.232E-02 

-0.275E-02 

-0.240E-02 

-0.173E-02 

-0.134E-02 

-0.47SE-03 


O.OQOE+OQ* 

O.OOOE+OO* 

Q.QQQE+00* 

O.OOOE+OO* 

O.OOOE+OO* 

O.OOOE+OO* 

O.OOOE+OO* 

O.OOOE+OO* 

O.OOOE+OO* 

O.OOOE+OO* 

O.OOOE+OO* 


0 -1 10 
ID 2 1/ 1/ 1 
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12 

0.170E-01 -0.133E-02 

O.OOOE+OO* 

13 

0.186E-01 -0.173E-02 

O.OOOE+OO* 

14 

0.170E-01 -0.214E-02 

0.000E+00* 

15 

0.154E-01 -0.299E-02 

Q.OOCE+QO* 

16 

0.136E-01 -0.213E-02 

O.OOOE+OO* 

17 

0.122E-01 -0.173E-02 

O.OOOE+OO* 

18 

0.130E-01 -0.110E-02 

O.OOOE+OO* 

13 

0.154E-01 -0.584E-03 

O.OOOE+OO* 

20 

0.177E-01 -0.109E-02 

O.OOOE+OO* 

21 

0.185E-01 -0.173E-02 

O.OOOE+OO* 

22 

0.177E-01 -0.238E-02 

O.OOOE+OO* 

23 

0.154E-01 -0.288E-02 

O.OOOE+OO* 

24 

0.130E-01 -0.237E-02 

O.OOOE+OO* 

25 

0.307E-01 -0.348E-02 

O.OOOE+OO* 

26 

0.264E-01 -0.340E-02 

O.OOOE+OO* 

27 

0.219E-01 -0.340E-02 

O.OOOE+OO* 

28 

0.174E-01 -0.371E-02 

O.OOOE+OO* 

29 

0.154E-01 -0.379E-02 

O.OOOE+OO* 

30 

0.134E-01 — 0 .371E-02 

O.OOOE+OO* 

31 

0.882E-02 -0.340E-02 

O.OOOE+OO* 

32 

0.430E-02 -0.341E-02 

O.OOOE+OO* 

33 

O.OOOE+OO* -0.347E-02 

O.OOOE+OO* 

34 

Q.307E-01 -0.311E-02 

O.OOOE+OO* 

35 

0.264E-01 -0.303E-02 

O.OOOE+OO* 

36 

0.220E-01 -0.300E-02 

O.OOOE+OO* 

37 

0.174E-01 -0.332E-02 

O.OOOE+OO* 

38 

0.154E-01 — 0.343E-02 

O.OOOE+OO* 

39 

0.133E-01 -0.332E-02 

O.OOOE+OO* 

40 

0.873E-02 -0.300E-02 

O.OOOE+OO* 

41 

0.433E-02 -0.303E-02 

O.OOOE+OO* 

42 

O.OOOE+OO* -0.310E-02 

O.OOOE+OO* 

43 

0.307E-01 -0.233E-02 

O.OOOE+OO* 

44 

0.264E-01 -0.232E-02 

O.OOOE+OO* 

45 

0.222E-01 -0.226E-02 

O.OOOE+OO* 

46 

0.853E-02 -0.226E-02 

O.OOOE+OO* 

47 

0.435E-02 -0.232E-02 

O.OOOE+OO* 

48 

0.000E+00* -0.233E-02 

O.OOOE+OO* 

49 

0.307E-01 -0.173E-02 

O.OOOE+OO* 

50 

0.264E-01 -0.173E-02 

O.OOOE+OO* 

51 

0.221E-01 -0.173E-02 

O.OOOE+OO* 

52 

0.867E-02 — 0 .173E-02 

O.OOOE+OO* 

53 

0.430E-02 -0.173E-02 

O.OOOE+OO* 

54 

0.000E+00* — 0 .173E-02 

O.OOOE+OO* 

55 

0.307E-01 -0.113E-02 

O.OOOE+OO* 

56 

0.264E-01 -0.115E-02 

O.OOOE+OO* 

c-? 

wi 

U.222E-01 -0.121E-02 

O.OOOE+OO* 

58 

0.353E-02 -0.121E-02 

O.OOOE+OO* 

59 

0.435E-02 -Q.115E-02 

O.OOOE+OO* 

60 

0.000E+00* -0.114E-02 

O.OOOE+OO* 

61 

0.307E-01 -0.362E-03 

O.OOOE+OO* 

62 

0.264E-01 -0.438E-03 

O.OOOE+OO* 

63 

0.220E-01 -0.469E-03 

O.OOOE+OO* 

64 

0.174E-01 -0.151E-03 

O.OOOE+OO* 
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65 

0.154E-01 -0.341E-04 

Q.OQQEiOQ* 

66 

0.133E-Q1 -0.153E-03 

O.OOOEiOO* 

67 

0.373E-O2 -0.469E-O3 

O.OOOEiOO* 

68 

O.433E-02 -0.433E-O3 

O.OOOEiOO* 

69 

0.0OOE+OO* -0.372E-O3 

O.OOOEiOO* 

70 

O.307E-O1 

0.125E-04 

O.OOOEiOO* 

71 

0.264E-01 -O.640E-O4 

O.QQOEiGQ* 

72 

0.213E-O1 -O.701E-O4 

O.OOOEiOO* 

73 

0.174E-O1 

Q.243E-03 

O.OOOEiOO* 

74 

O.154E-01 

0.325E-03 

O.OOOEiOO* 

75 

0.134E-01 

0.242E-03 

O.OOOEiOO* 

76 

O.882E-02 -0.701E-04 

O.OOOEiOO* 

77 

0.430E-02 -0.596E-O4 

O.OOOEiOO* 

78 

O.OOOEiOO* O.OOOEiOO* 

O.OOOEiOO* 


OLIB READ ERROR, NU,l,KORE,IERR,IOP= 1 
NAME* CASE TITL 11263747405 

1 STATIC REACTIONS, FORCE ERRORS. 


OJQINT 

1 

2 

3 

1 

0.100E-06 

0.109E-06 

O.OOOEiOO* 

2 

0.583E-06 -0.224E-06 

O.OOOEiOO* 

3 

0.982E-06 

0.220E-06 

O.OOOEiOO* 

4 

0.659E-06 -0.384E-06 

O.OOOEiOO* 

5 

-0.584E-06 

0.113E-05 

O.OOOEiOO* 

6 

0.177E-06 

0.668E-06 

O.OOOEiOO* 

7 

0.695E-06 -0.108E-06 

O.OOOEiOO* 

8 

0.352E-06 -0.740E-06 

O.OOOEiOO* 

9 

0.140E-05 -0.858E-Q7 

O.OOOEiOO* 

10 

-0.967E-06 -0.569E-06 

O.OOOEiOO* 

11 

-0.231E-05 -0.168E-06 

O.OOOEiOO* 

12 

-0.252E-05 -0.122E-05 

O.OOOEiOO* 

13 

-0.177E-05 

0.186E-05 

O.OOOEiOO* 

14 

-0.456E-05 -0.245E-08 

O.OOOEiOO* 

15 

-0.491E-06 -0.823E-06 

O.OOOEiOO* 

16 

0.240E-05 

0.284E-06 

O.OOOEiOO* 

17 

-0.113E-05 -0.833E-06 

O.OOOEiOO* 

18 

0.318E-05 -0.154E-06 

O.OOOEiOO* 

19 

0.867E-06 

0.162E-06 

O.OOOEiOO* 

20 

0.411E-05 -0.230E-06 

O.OOOEiOO* 

21 

0.488E-05 -0.138E-05 

O.OOOEiOO* 

22 

0.213E-05 

0.139E-05 

O.OOOEiOO* 

23 

0.370E-Q6 -0.750E-06 

O.OOOEiOO* 

24 

-0.323E-06 

0.103E-05 

O.OOOEiOO* 

25 

0.954E-06 

0.149E-06 

O.OOOEiOO* 

26 

0.137E-06 -0.454E-06 

O.OOOEiOO* 

27 

0.974E-06 

0.388E-06 

O.OOOEiOO* 

28 

0.414E-06 

0.415E-06 

O.OOOEiOO* 

29 

-0.474E-06 

0.488E-06 

O.OOOEiOO* 

30 

0.525E-06 

0.123E-06 

O.OOOEiOO* 

31 

-0.307E-06 

0.436E-06 

O.OOOEiOO* 

32 

-0.538E-07 

0.813E-07 

O.OOOEiOO* 

33 

-0.620Ei00* -0.365E-07 

O.OOOEiOO* 

34 

-Q.238E-Q6 

0.453E-06 

O.OOOEiOO* 

35 

0.128E-05 -0.149E-06 

O.OOOEiOO* 

36 

0.350E-05 -0.367E-06 

O.OOOEiOO* 


37 0.240E-05 -0.118E-05 O.OOOE+OO* 

38 -Q.346E-05 -0.531E-06 O.OOOE+OO* 

39 0.154E-05 -0.723E-07 O.OOOE+OO* 

40 -Q.751E-06 0.330E-06 O.OOOE+OO* 

41 -0.959E-07 -0.841E-06 O.OOOE+OO* 

42 -0.137E+01* -0.146E-G6 O.OOOE+OO* 

43 -0.238E-05 0.464E-06 O.OOOE+OO* 

44 0.275E-05 0.106E-05 O.OOOE+OO* 

45 -O.655E-05 -0.226E-06 O.OOOE+OO* 

46 0.746E-06 -0.305E-06 O.OOOE+OO* 

47 Q.911E-06 0.415E-06 O.OOOE+OO* 

43 -0.225E+01* 0.299E-06 O.OOOE+OO* 

49 -0.238E-05 0.104E-05 O.OOOE+OO* 

50 -0.330E-05 -0.676E-Q6 O.OOOE+OO* 

51 -0.697E-06 -0.187E-Q5 O.OOOE+OO* 

52 -0.477E-06 0.201E-06 O.OOOE+OO* 

53 -0.10QE-05 0.341E-06 O.OOOE+OO* 

54 -0.202E+01* 0.592E-07 O.OOOE+OO* 

55 0.334E-05 -0.133E-06 O.OOOE+OO* 

56 -0.287E-05 0.136E-Q5 O.OOOE+OO* 

57 -0.159E-05 0.958E-06 O.OOOE+OO* 

58 -0.104E-05 -0.367E-06 O.OOOE+OO* 

59 -0.202E-06 -0.314E-06 O.OOOE+OO* 

60 -0.225E+01* -0.458E-07 O.OOOE+OO* 

61 0.715E-06 -0.646E-06 O.OOOE+OO* 

62 0.193E-05 -Q.162E-05 O.OOOE+OO* 

63 Q.291E-05 0.101E-05 O.OOOE+OO* 

64 -0.752E-06 0.349E-06 O.OOOE+OO* 

65 -0.190E-07 0.553E-06 O.OOOE+OO* 

66 0.436E-06 -Q.145E-06 O.OOOE+OO* 

67 0.716E-07 -0.438E-06 O.OOOE+OO* 

68 -0.166E-06 0.946E-07 O.OOOE+OO* 

69 -0.187E+01* -Q.121E-06 O.OOOE+OO* 

70 -Q.119E-06 0.276E-07 O.OOOE+OO* 

71 -0.135E-05 0.967E-06 - O.OOOE+OO* 

72 -0.175E-05 0.204E-06 O.OOOE+OO* 

73 0.854E-06 -0.544E-06 O.OOOE+OO* 

74 -0.754E-06 -0.336E-07 O.OOOE+OO* 

75 -0.101E-Q5 -0.187E-06 O.OOOE+OO* 

76 -0.411E-06 -0.285E-06 O.OOOE+OO* 

77 -0.190E-06 -0.155E-06 O.OOOE+OO* 

78 -0.620E+00* 0.493E-05* O.OOOE+OO* 

EXIT VPRT CPUTIHE= 3.3 I/0(D!R,BUF)= 0 0 

<W1> CLOSE, Idi! 1, File; PAN1.L01 
<CL> CSS exhausted 

ENDRUN called by CLIP 
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APPENDIX C 


PROGRAM FOR ELASTO-PLASTIC STABILITY PROBLEM 

The program written in PASCAL is divided into tvx> parts. Part I is 

called PROGRAM COLUMN. This part of the program is down-loaded onto one 

of the NPR number of processors and this processor is designated as 
"Master Processor". Part II of the program is named PROGRAM TANGENT. 
This program is down-loaded onto the rest of the processors of the NPR 
number of processors. These (NPR-1) number of processors are called 

"Assistant Processors". The assistant processors are utilized to compute 

elastic-plastic cross-sectinal properties of various cross sections along 
the length of the member which are communicated to the master processor. 
The master processor assembles the global stiffness matrix upon receipt of 
the cross sectional properties form assistants. The lateral displacement 
vector is computed on the master processor. Thus the operation is 
sequential between the assembly of global stiffness values and evaluation 
of the lateral displacements of the member. A detailed flow chart of the 
algorithm used and the concurrent computational procedure is outlined in 
Reference 6. 

The input to the program consists of the member length, material 
constants, end rotational stiffness characteristics, initial out-of- 
straightness midspan amplitudes, cross-sectional dimensions, and residual 
stresses. The output from the program consists of axial load level, 
corresponding converged lateral displacements of the member, 

nondimensionalized determinant of the global stiffness matrix at each load 
level, and the computational times of each processor. 
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ORIGINAL PAGE JS 
OF POOR quality: 

PROGRAM FOR ELASTO-PLASTIC STABILITY PROBLEM 


< s- * W I DEL I ST , NO ASSERTS , NO TRACEBACK * ) 

<* 3GL0BAL0PT *> 

PROGRAM COLUMN; 

CONST MAX I DX =255? <*MAX NBR I/O INDEX TAG *) 

M A X REC=255 * <*MAX REC LENGTH FOR NBR I/O *) 
MAXN0DE=36? <*MAX NODE NUMBERS *) 

MAXDA=31? ( -s-MAX DATA AREA NUMBER #> 

MAX I NT =32767 * <*MAX INTEGER -a-) 

SYSFLAG-1 ? (^SYSTEM FLAG *> 

TYPE NODE® 1. . . MAX NODE? 

I D X = 1 . . MAXIDX? 

RECLEN- 1 . . MAX REC? 

DANUM= 1 . . MAXDA? 

FLAG=0. . 7; 

ADDR= INTEGER; 

F’OS I NT — 1 . .MAX I NT? 

<* PASLIB EXTERNAL PROCEDURES AND FUNCTIONS FOLLOW *) 


(* TIMER PROCEDURES — *) 

FUNCTION XTIME? LONG I NT? EXTERNAL; 

PROCEDURE T ST ART ( T s POS I NT ) ? EXTERNAL ? 

FUNCT I ON TREAD : LONG I NT ? EX TERNAL ? 

PROCEDURE TSTOP ? EXTERNAL ? 

<* TEXT OUTPUT PROCEDURES *> 

PROCEDURE MSG ( STR I MG : PACKED ARRAY Cl. . ? I OF CHAR ) ? E X TERNAL ? 
PROCEDURE MSGL.N ( STR I NO s PACKED ARRAY Ct. . ? 1 OF CHAR ) ? E X i ERNAL 
PROCEDURE NXTLN? EXTERNAL? 

PROCEDURE MSG I ( I : INTEGER) ? EXTERNAL? 

PROCEDURE MSGR ( X s REAL ) ? EXTERNAL ? 

PROCEDURE MSGL ( I : LONG I NT ) ? EXTERNAL ? 

PROCEDURE ENDLN ( N : POS I NT ) ? EXTERNAL ? 

( ?<• >- } 

(•» INTERACTIVE INPUT *) 

( * Vr ) 

PROCEDURE QUERY ? EXTERNAL ? 

F! JNCT I ON RD I s I NTEOER ? E X TERNAL ? 

FI JNCT I ON RDR : REAL ? EXTERNAL ? 

( s # ) 

FUNCTION DAPTR ( DA ? DANUM) : ADDR? EXTERNAL? 
f <- a- ) 

( * FLOATING OPERATIONS * ) 

FI IMC T T ON ADD ( X , Y : REAL ) : REAL ? EXTERNAL ? 
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FUNCTION 

FUNCTION 

FUNCTION 

FUNCTION 

FUNCTION 

FUNCTION 

FUNCTION 

<* 

< * 

FUNCTION 

FUNCTION 


SUB ( X , V : REAL ) : REAL; EXTERNAL? 
MULT < X , Y : RE AL ) : REAL ; EXTERNAL ; 
DI VD ( X , Y: REAL ) : REAL; EXTERNAL; 
NEG ( X : REAL ): REAL ; EXTERNAL ; 
ABS95 ( X : REAL ) : REAL; EXTERNAL ; 
CMP < X , Y : REAL ) : REAL ; EXTERNAL ; 
SQRT95 < X : REAL ) ; REAL ; EXTERNAL ; 


— FLOATING POINT CONSTANTS 
MAX95: REAL; EXTERNAL; 

MIN95! REAL; EXTERNAL; 


*> 

*) 


<* FLOATING POINT CONVERSIONS *) 

FUNCTION IFI X ( X s REAL ) ! INTEGER; EXTERNAL; 

FUNCT ION FLO AT I ( I : I NTEGER ) : RE AL ! E X TERN AL ; 


I i 1.1 1 M >_ 


FUNCTION VDP(N:POSINT; VAR A: ARRAYC 1 . . ?] OF REAL; 

VAR B: ARRAY C 1 . . ?] OF REAL) : REAL; EXTERNAL; 

<* NBR COMMUNICATIONS *> 

PROCEDURE SEND2(N: NODE; INDEX; IDX; LOC: ADDR; NWORDSs RECLEN) ; EXTERNAL; 
PROCEDURE SEND2ALL( INDEX: I DX; LOC; ADDR; NWORDS: RECLEN ); EXTERNAL; 

procedure recv2(n; node; index; idx; loc;addr; nwords; reclen); external; 

<* FLAG ROUTINES *) 


procedure flgen ( f; flag > ; external; 

PROCEDURE FLGSET (F: FLAG ); EXTERNAL; 
PROCEDURE FLORES (F; FLAG) ; EXTERNAL; 


PROCEDURE BAR ( F : FLAG ) ; EXTERNAL ; 

FUNCTION ANY (F: FLAG ) : BOOLEAN; EXTERNAL; 

<* PROCESSOR IDENTIFICATION *> 


PROCEDURE PSCL$$; 

CONST 

EDIM=3; 

TYPE 

DA7= ARRAY Cl. .33 OF I NTEGER ; 

DA8=ARRAYC1. .4] OF REAL; 

VAR 

pb;real; nentry; integer; data; array ci..4] of real; 

NEW, IO,NPR,NSEC,NPROB,N,MDIM ; INTEGER; 

ELASTIC : BOOLEAN; 

Z ERO , ONE , TWO : RE AL ; 

DET I , DET , P I NC , TOL 1 , T0L2 , TOL3 : REAL ? 

DALL7 s@DA7; DALL8 : <5DA3; 

CLK1 » CLK3 : LONG I NT; 

PROCEDURE GLOBAL ( N , MD I M , NPRGB , NSEC : I NTEGER ) ; 
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LABEL 1 0 i 20 


TYPE 

VECT0R=ARRAYC 1 . . N3 OF REAL? 

MATRI X=ARRAYC l . . N, 1 . . N3 OF REAL? 

VAR 

K:MATRIX; D I S , D I » D , F : VECTOR ; 

UV: ARRAYC1. .MDIM3 OF REAL? 

COEF : ARRAY Cl. . 143 OF REAL? IDSEC, ITER, TANIDX: INTEGER? 

COUNT, CNTR, NITER, I SEC: INTEGER; 

READY: ARRAY C 1 . . NPR, 1 . . NPR0B3 OF BOOLEAN; 

I CNT : INTEGER; 

( tt*##*###*#########***#*##**#*#**###*##-*#*-***#*#-*#***##*#*## ) 

PROCEDURE MATINV< VAR A: MATRIX; N: INTEGER; VAR DET: REAL ) ; 

<* MATRIX INVERSION PROCEDURE *> 

(* METHOD - JOHNSON ( 1 ) METHOD #) 

<* A IS N X N SYMMETRIC MATRIX#) 

VAR 

IPV: ARRAYC1. .N, 1. .33 OF INTEGER; 

NSW AP , L , I , I COL , JCOL, I ROW , JROW : I NTEGER ; 

TEMP, T,PIV, SWAP, AMAX : REAL; 


BEGIN 

DET:=ONE; TEMP: =AC 1 , ID; 

FOR J:=l TO N DO IPVCJ,33:=0; 

I : =1 ; 

REPEAT (* UNTIL DET=0 OR I>N *) 

AMAX : =ZERQ; 

FOR J: =1 TO N DO 
BEGIN 

IF IPVC J, 3101 THEN 
BEGIN 

FOR KK : = 1 TO N DO 
BEGIN 

IF IPVCKK, 3301 THEN 

IF CMP < AMAX , ABS95 ( AC J, KK3 ) ) <0 THEN 
BEGIN 

I ROW : = J ; <* SAVE ROW & COLUMN *> 
ICOL:=KK; 

AMAX : =ABS95 < A C J , KK 3 ) ; 

END; <* IF AMAX #) 

END; (* FOR KK *> 

END; <* IF IPVC J, 33 *) 
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END; <*• FOR J *) 

(* IF PIVOT ELEMENT IS EQ 0 - DET = 0 *> 
IF CMP ( AMAX » T0L2 ) <0 THEN DET:=ZERO 


ELSE 

BEGIN 

IPVC I COL, 33 s =IPVC I COL, 33+1 * 

ipvci, i]:=irow; 
ipvc I » 23 ; =icol; 

<* INTERCHANGE ROWS TO PUT PIVOT ELEMENT ON DIAGONAL *> 

IF IROWOICOL THEN 

FOR j;»l TO N DO 
BEGIN 

SWAP: =A Cl ROW, J3; 

AC IROW, J 3 5 = A C I COL, J3 ; 

A C I COL , J 3 s =SW AP ; 

END; <» FOR J 1 TO N *> 

<* DIVIDE PIVOT ROW BY PIVOT ELEMENT *> 

PIV: =AC I COL, IC0L3; 

DET : =MULT ( DET , D I VD < P I V , TEMP ) ) ; 

ACICOL, IC0L3S=0NE; 

FOR KK : = 1 TO N DO 

A C I COL , KK 3 : =D I VD ( A C I COL , KK 3 , P I V ) ; 

<* REDUCE THE NON-PIVOT' ROW BY SUBTRACTION *> 

FOR Ll:«l TO N DO 

IF LIOICOL THEN (* AVOID DIAGONAL *> 

BEGIN 

T : =A C L 1 , I COL 3 ; 

ACL1 , IC0L3 s =ZERO» 

FOR KK*.=1 TO N DO 

A C L 1 , KK 3 : =SUB ( A C L 1 , KK 3 , MULT ( A C I COL , KK 3 , T ) ) 

END; <* IF LI O I COL *) 

I ; =1+1 ; <* GO ON TO NEXT *) 

END; <* ELSE IF AMAX-CTOL *> 

UNTIL ( I >N ) OR CMP (DET, ZERO )*0» 

<* INTERCHANGE THE COLUMNS AND MODIFY DET *> 

NSWAPs =0; <* THIS IS SIGN FLAG *> 

IF CMP (DET, ZERO) <>0 THEN BEGIN (* DO ONLY IF DET. NE. 0.0 #) 
FOR J: =1 TO N DO 
BEGIN 
L:=N-.J+l; 

IF I P V C L , 13 <> I P V C L , 2 3 THEN 
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BEGIN 

JROW:=IPVCL, 13? 

JCOL s = I P V C L , 2 3 ; 

NSWAP: =NSWAP+1 ; (* COUNT SWAPS *) 

FOR KK ! = 1 TO N DO 
BEGIN 

SWAP : = A C KK , JROW 3 ; 

A C KK , JROW 3 : = A C KK , JCOL 3 ? 

A C KK , JCOL 3 s =SW AP * 

END? <* FOR KK *> 

END; <* IF IPV *) 

END; (* FOR J *) 

IF ODD (NSWAP) THEN SWAP: =NEG (ONE) 

ELSE SWAP: =ONE; 

END; (* IF DETOO.O *) 

DET : =MULT ( DET , SWAP ) ; 

END; <* MATINV PROCEDURE *) 

( ) 
PROCEDURE SOLVE ( N: INTEGER; VAR X: MATRIX; VAR SOLNS: VECTOR; 

SUMS: VECTOR; VAR DET: REAL); 

BEGIN (» MAIN PROGRAM *> 

MATINV (X,N, DET); (* CALL MATINV PROCEDURE *) 

IF CMP (DET, ZERO) =0 THEN 

MSGLN< ' MATRIX IS SINGULAR' > 

ELSE 

BEGIN 

FOR J: =1 TO N DO 
BEGIN 

SOLNS C J 3 : =ZERO ; (* INITIALIZE *) 

SOLNS C J 3 : = VDP ( N , X C J 3 , SUMS ) ; 

END; <* FOR J *) 

END; (* ELSE i.«. DET .NE. 0.0 #) 

END; 

PROCEDURE ASSEMBLE ( I DSEC : I NTEGER ) ; 

VAR IK: INTEGER; 

BEGIN 

(*$ NO OPTIMIZE *) 

IF IDSEC03 THEN 
BEGIN 

IF I DSEC= 1 THEN 
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BEGIN 

K C l , 13 : = ADD ( COEF C 1 3 , COEF C 5 3 ) ! 

K C 1 > 2 3 s =COEF C 3 3 ' 

K C 1 , 3 3 : =SUB < COEF C 1 3 , COEF C 5 3 ) ? 

K C 1 7 4 3 s =COEF C 3 3 ? 

KC 1 , N-l 3 : =NEG ( COEF C 63 ) ! 

K C 1 , N-3 3 ! =COEF C 6 3 * 

FC 1 3 5 =COEFC73 ? 

KC2, 1 3 ! =COEF C 8 3 ; 

K C 2 , 2 3 : = ADD ( COEF C 1 0 3 , COEF C 123 ) 5 
K C 2 1 3 3 s =COEF C S 3 * 

K C 2 , 4 3 : =SUB ( COEF C 1 0 3 , COEF C 1 2 3 ) 5 
K C 2 7 N 3 : =NEG ( COEF C 1 3 3 ) ! 

K C 2 7 N-2 3 : =COEF C 1 3 3 ? 

F C 2 3 ! =COEF C 1 4 3 * 

END 

EL3E 

IF IDSEC=2 THEN 
BEGIN 

K C 3 7 13: =COEF C 5 3 ; 

K C 3 , 3 3 : -SUB ( COEF C 2 3 , COEF C 5 3 ) ; 
KC3> 43 ! =COEFC43 * 

K C 3 7 5 3 s =COEF C 1 3 * 

K C 3 7 6 3 : =COEF C 6 3 ! 

K C 3 7 N- 1 3 s =NEG < COEF C 6 3 ) ; 

K C 3 7 N-3 3 ! =COEF C 6 3 * 

F C 3 3 s =COEF C 7 3 ? 

K C 4 7 2 3 ! =COEF C 1 2 3 * 

K C 4 , 3 3 ! =COEF C 9 3 * 

K [ 4 , 4 3 : =SIJB ( COEF C 1 1 3 , COEF C 123 > 5 
K C 4 , 5 3 : =COEF C 8 3 * 

K C 4 , 6 3 ! =COEF C 1 0 3 ! 

K C 4 7 N 3 : =NEG ( COEF C 1 3 3 ) ? 

K C 4 , N-2 3 J =COEF C 1 3 3 * 

F C 4 3 s =COEF C 1 4 3 * 

END 

ELSE 

BEGIN 

K C 5 7 13s =COEF C 5 3 5 
K C 5 7 3 3 : =SUB ( COEF C 1 3 7 COEF C 5 3 ) ! 
K C 5 , 4 3 5 =COEF C 3 3 ? 

K C 5 , 5 3 : =COEF C 2 3 ; 

K C 5 7 6 3 : =COEF C 4 3 ! 

K C 5 , 7 3 s =COEF C 1 3 * 

K C 5 7 8 3 : =COEF C 3 3 5 



K E 5 , N- 13 : =NEG ( COEF E 6 3 ) ; 

KE5, N-33 ! =CQEFE63 ? 

F E 5 3 ! =COEF E 7 3 ; 

KE6, 23 s =COEFE 123 * 

K E 6 , 3 3 s =COEF E 8 3 » 

K E 6 , 4 3 : -SUB ( COEF E 1 0 3 , COEF E 1 2 3 ) ; 

K E 6 , 5 3 : =COEF E 9 3 > 

K E 6 » 6 3 s =COEF E 1 1 3 * 

K E 6 , 7 3 s =COEF C 8 3 * 

K E 6 » 8 3 s =COEF E 1 0 3 * 

K C 6 1 N 3 : =NEG ( COEF E 1 3 3 ) ; 

KE6, N— 2 3 : -COEF E 1 3 3 ? 

FC63 s =COEFE 143 ; 

END? <* I DSEC=3 *> 

END <* I DSEC03 THEN *> 

ELSE 

IF IDSEO-NSEC-2 THEN 
BEGIN 

IKs = < EDIM—1 ) *IDSEC; 

IF I DSEC=NSEC— 2 THEN 
BEGIN 

KE IK-1, 13: =COEFE53; 

KE IK-1, 33 :=NEG(COEFE53); 

KE IK-1 , N— 73 ; =COEF E 1 3 ; 

K C I K- 1 , N-6 3 : =COEF E 3 3 ; 

K E I K- 1 , N— 5 3 : =COEF C 2 3 ; 

K E I K- 1 , N— 4 3 : =COEF C 4 3 ; 

K E I K- 1 , N-3 3 : = ADD ( COEF C 1 3 , COEF E 6 3 ) ; 
K C I K- 1 , N-2 3 : =CQEF C 3 3 ; 

K C I K- 1 , N- 1 3 : =NEG ( COEF E 6 3 ) ; 

FE IK-1 3 ; =COEFE73 * 

KE IK, N— 73 s =COEFE83 * 

K E I K , N-6 3 : =COEF E 1 0 3 ; 

KE IK, N— 53 s =COEFE93 ; 

KE IK , N— 4 3 : =COEF E 1 1 3 ; 

KE IK, N-33 ; =COEFE83 * 

K E I K , N— 2 3 : = ADD ( COEF E 1 0 3 , COEF E 1 3 3 ) ; 
K E I K , N 3 s =NEG ( COEF E 1 3 3 ) ? 

KE IK, 23 s =COEFE 123 * 

KE IK, 43 : =NEG ( COEF E 123 ) ? 

F C I K 3 : =COEF E 1 4 3 * 

END 
ELSE 

IF IDSEC-NSEC- 
BEGIN 


1 THEN 


K C IK-1 , N— 5 3 : =COEF E13; 

K Z I K- 1 , N-4 3 : =CQEF C 3 3 ! 

K C I K- 1 , N-3 3 s = ADD ( COEF C 2 3 , COEF C 6 3 ) ? 

K C I K- 1 , N-2 3 : =CQEF Z 4 3 ? 

KC IK-1 , N-l 3 ! =NEG< COEFC63 ) ? 

KC IK-1, U:=COEFC53; 

KC IK-1, 33: =NEG ( COEF C 5 3 ) ; 

FC IK-13 s =COEFC73 * 

KC IK, N-53 s =COEF C 8 3 ; 

K C I K , N— 4 3 * =COEF C 1 0 3 * 

K C I K , N-3 3 : =COEF C 9 3 ; 

K C I K , N-2 3 : = ADD ( COEF C 1 1 3 , COEF C 1 3 3 ) U 
KC IK, N3 : =NEG ( COEF L 1 3 3 ) 5 
KC IK, 23 s =COEFC 123 5 
K C I K , 4 3 : =NEG ( COEF C 1 2 3 ) ? 

FC IK3 s =COEFC 143 ? 

END 

ELSE 

BEGIN 

KC IK-1 , N-3 3 s « ADD (COEF C 1 3 , COEF C 63 ) ; 
K C I K- 1 , N-2 3 : =COEF C 3 3 ? 

K C I K- 1 , N- 1 3 : =SUB ( COEF C 1 3 , COEF C 6 3 ) ? 
KC IK-1 , N3 s =COEFC33 5 
KC IK-1, 13s =COEF C 5 3 * 

KC IK-1 , 33 s =NEG( COEFC53 > ! 

FC IK-1 3 s =COEFC73 * 

K C I K , N-3 3 : =COEF C 8 3 ? 

K C I K , N-2 3 : = ADD ( COEF C 1 0 3 , COEF C 1 3 3 ) ! 
KC IK, N— 1 3 s =COEFC83 * 

K C I K , N 3 : -SUB ( COEF C 1 0 3 , COEF C 1 3 3 ) 5 
KC IK, 23 ! =COEFC 123 * 

K C I K , 4 3 : =NEG ( COEF C 123 ) 5 
F C I K 3 ! =COEF C 1 4 3 * 

END? 

END <* I DSEC>=NSEC-2 THEN *) 

ELSE 

BEGIN 

IK! =<EDIM-1 ) #IDSEC? 

KC IK-1, 13!=COEFC53! 

K C I K- 1 , 3 3 : =NEG ( COEF C 5 3 > ? 

K C I K- 1 , N-3 3 : =COEF C 6 3 ? 

K C I K- 1 , N- 1 3 ! =NEG ( COEF C 6 3 ) 5 
K C I K- 1 , 2* ( I DSEC- 1 ) - 1 3 ! =COEF C 1 3 ? 

K C IK-1 , 2* < I DSEC- 1)3! =COEF [33? 

K C I K- 1 , 2* I DSEC- 1 3 ! =COEF C 2 3 ! 


K C I K- 1 , 2* I DSEC 3 : =COEF C43; 

K C I K- 1 , 2* ( I DSEC+1 ) - 1 3 : -COEF 111 ’, 

K C I K- 1 , 2* ( I DSEC+ 1)3: -COEF C 3 3 ; 

F C I K- 1 3 : —COEF C 7 3 ; 

KC IK, 23 : -COEFC 12] * 

K C I K , 4 3 s -NEG ( COEF C 1 2 3 ) " 

K C I K , N-2 3 : -COEF C 1 3 3 ; 

K C I K i N 3 : =NEG ( COEF C 1 3 3 ) ! 

K C I K j 2* < I DSEC- 1 ) - 1 3 : =COEF C 8 3 5 
KC IK, 2* < I DSEC- 1)3: -COEF C 1 0 3 5 
KC IK, 2*1 DSEC- 1 3 : -COEF C 9 3 ; 

KC IK, 2*1 DSEC 3 : -COEF C 1 1 3 S 
K C I K , 2* < I DSEC+ 1 ) - 1 3 : -COEF C S 3 5 
KCIK,2*( I DSEC+1 ) 3:=C0EFC103; 

F C I K 3 : =COEF C 1 4 3 * 

END; (# I DSEC other than (<=3) O-NSEC-2) * ) 

END; (* Assemble end #) 

BEGIN 

NENTRY: =0; ICNT:=0; CNTR:=0; CLK3:=0; 

ELASTIC: -TRUE; 

FOR I : = 1 TO N DO 
BEGIN 

FOR J : - ' 1 TO N DO KC I , J3 : -ZERO; 

DC I 3 : -ZERO; 

DI C I 3 : -ZERO; 

FC I 3 : -ZERO; 

KC I , 13: -ONE 
END; 

FOR I : =1 TO MDIM DO UVCI3:=ZER0; 

REPEAT SMALL PINC *) 

REPEAT (* UNTIL CAPACITY *> 

PB:=ADD(PB,PINC>; 

NITER: =0; 

REPEAT (* UNTIL UV CONVERGES *) 

DAT A C 1 3 : -SUB (DC43,DC23>; 

DATA C 2 3 : -SUB ( D C N-2 3 , D C N 3 ) ; 

DATA C 3 3 : -SUB ( D C 3 3 , D C 1 3 ) ; 

DAT A C 4 3 : -SUB ( D C N-3 3 , D C N- 13); <* Terms to set Mbx, Mtx, . . etc . # ) 

CLK 1 : -TREAD ; SEND2 ALL ( 1 , LOCAT I ON ( PB ) , 1 1 ) ; CLK3 : -CLK3+TRE AD-CLK 1 ; 
NEW: =1 ; 

I SEC: =0; 

FOR I : = 1 TO NPR DO 
BEGIN 

FOR J: - 1 TO MDIM DO UVC J3 : -DC ISEC+J3 ; 
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I SEC s = I SEC+MD I M * CLK 1 ? =TREAD? 

SEND2 <1,2, LOCAT I ON ( UV ) , 2*MD I M ) ? CLK3 ? -CLK3+TRE AD-CLK 1 
END? 

CLK1 : =TREAD? 

BAR < 2 ) ? 

CLK 3 ? =CLK3+TRE AD-CLK 1 ? 

FOR I : =1 TO N DO FOR J: = l TO N DO KC I » JH ? =ZERO? 

FOR is *1 TO NPR DO FOR J?=l TO NPROB DO READY C I , J] : =FALSE? 
COUNT? =0? 

WHILE (COUNT<NSEC> DO <* Asynchronous Communication *) 
FOR I?= 1 TO NPR DO 
FOR J? = 1 TO NPROB DO 
IF NOT READY C I » J ] THEN 
BEGIN 

i~ LK 1 5 =TREAP? 

RECV2 < I , J, LOCATION < COEF) , 31 ) ? CLK3? -CLK3+TRE AD-CLK 1 ? 

IF ITERMCNT THEN 
BEGIN 

IF TANIDX=i THEN GOTO 20? 

READY C I , J 3 ? =TRUE ? 

COUNT ?=C0UNT+1? 

IF I0<=5 THEN MSOKIDSEC); 

ASSEMBLE < I DSEC) 

END? 

END? 

IF I0<=3 AND ANY (5) THEN 
BEGIN 

NXTLN? 

FOR I : =1 TO N DO BEGIN 

FOR J ? = 1 TO N DO BEG I N MSGR < K L I , J 1 ) ? MSG < ' ' ) END ? 

ENDLN < 2 ) ? END? 

END? 

IF NENTRY=0 THEN 
BEGIN 

SOLVE <N,K,D,F,DET>? 

DETI ? =DET? DET? =DIVD(DET, DETI ) ? 

NENTRY? =1 ? 

FOR is—l TO N DO DICn?=FCI]? 

END 

ELSE 

BEGIN 

SOLVE <N,K,D,F, DET ) ? 

DET? =DIVD< DET, DETI)? 

IF 1004 THEN BEGIN NXTLN? MSGR < DET )? NXTLN? END? 

IF CMP < DET , T0L2 ) <=0 THEN GOTO 20? 
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CNTRs =0? 

FOR I : = 1 TO N DO 
BEGIN 

I F CMP ( ABS95 ( SUB ( F t 13 , D I C 1 3 ) ) » TOL 1 ) OO THEN CNTR : =CNTR+ 1 ? 

END 5 
END? 

IF (CNTRON) THEN FOR I: =1 TO N DO DICI3:=FCI3; 

IF NOT ANY (5) THEN NEW:=0 ELSE IF CNTR=N THEN NEW:=0 ELSE NEW:=l; 

CLK 1 ; =TRE AD ; SEND2 ALL ( 3 , LOCAT I ON ( NEW ) , 1 ) 5 BAR (7)5 CLK3 : =CLK3+TRE AD-CLK 1 ? 
ICNT:=ICNT+l; 

NITER: =NITER+1; 

UNTIL ( CNTR=N ) OR (NITER=8> 5 

IF NITER>=8 THEN BEGIN MSG ("NO GLOBAL CONVERGENCE " ) * GOTO 20 END; 

IF 1005 THEN 
BEGIN 

NX TLN ; MSG ("AXIAL LOAD= " ) ; MSGR ( PB ) ; 

NXTLN; MSGLN( "DISPLACEMENT VECTOR" ) ; 

FOR I : =1 TO N DO BEGIN MSGR ( DC 1 3 ) ; MSG ( " "); END; NXTLN; 

END; 

FOR I : =1 TO N DO DISC 1 3 : =DC I 3 ; 

IF I0>=7 THEN BEGIN MSG( "AXIAL LOAD=" >; MSGR < PB >; NXTLN; END; 

MSG( "Determinant of Global Stiffness Matrix =">; MSGR < DET ); NXTLN; 

UNTIL CMP(DET,T0L2)O0; 

20; • 

IF TANIDX = 1 OR CMP ( DET » T0L2 ) OO THEN ICNT: =ICNT+1 ; 

IF NITERC3 THEN 

IF NEW=1 THEN BEGIN CLK1 : =TREAD; 

SEND2ALL ( 3 , LOCAT I ON ( NEW ) , 1 ) ; BAR ( 7 ) ; 

CLK 3 : =CLK3+TREAD-CLK 1 ; END ; 

for i ; =1 to n do Dtn;=Discn; 

PB: =SUB(PB> PINO ; (* Retain non i ncremented PB *) 

P I NC : =D I VD ( P I NC , T WO ) ; 

UNT I L CMP ( P I NC , TOL 3 ) OO ; 

10: FLGSET ( 6 ) » 

BAR (2); 

NXTLN; 

END; (* PROCEDURE GLOBAL *) 


BEGIN (* MAIN PROGRAM *) 

FLGEN ( 2 ) ; FLGEN ( 5 ) ; FLGEN ( 6 ) ; FLGEN ( 7 ) ; 


FLGRES ( 2 ) ; FLORES ( 5 ) ; FLGRES (6)3 FLGRES ( 7 ) ; BAR ( SYSFLAG ) 


ZERO : =CV?5 12(0.0); ONE : =CV95 12(1.0); TWO : =CV95 12(2.0); 
DALL7 : : ADDR : =DAPTR ( 7 ) ; BALLS: : ADDR: =DAPTR ( S ) ; 

NSEC: =DALL7@C 1 3 ; NPR: =DALL7@C23 ; 10: =DALL7@C33 ; 

TOL 1 : =CV95 1 2 ( D ALLS® C 1 3 ) 5 T0L2 : =C V95 1 2 ( D ALLS© C 23); 
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T0L3 : =CV95 1 2 ( D ALL 80 C 3 1 > ! 

P T NO : =C V95 1 2 < BALLS® C 4 1 > 5 PB : = ZERO ; 

1ST ART ( 1 ) ; 

MSG ( •' CROSS SECT I ONS* ' ) 5 MSG I ( NSEC ) ; 

MSG ( - PROCESSORS* x > ! MSG I ( NPR+ 1 ) ; NXTLN 5 

NPROBs =NSEC DIV NPR5 
MB I M s =NPROB* < ED I M- 1 ) ; 

N:=NSEC*(EDIM-1>; 

GLOBAL (N,MDIM, NPROB , NSEC > 5 
TSTOP? 

ENDLN < 2) ; MSG < 'PROCESSING TIME *' ) ?MSGL( TREAD ) ; NXTLN; 
NXTLN? MSG ( ' EXECUTION TIME = ' ) ; MSGL ( TREAD-CLK3 ) ; NXTLN; 

end; 


BEGIN <**NO OBJECT*) 
END. 
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< *$WI DEL I ST , NO ASSERTS , NO TRACEBACK* ) 

( * $ GLOBAL OPT *) 

PROGRAM TANGENT? 

CONST MAX I DX =255? <*MAX NBR I/O INDEX TAG *> 

MAXREC=255* <*MAX REC LENGTH FOR NBR I/O *> 
MAXN0DE=36? <*MAX NODE NUMBERS ») 

MAXDA=3l! <*MAX DATA AREA NUMBER *> 

MAX INT=32767? (#MAX INTEGER *> 

SYSFLAG=l! (^SYSTEM FLAG *> 

TYPE NODE= 1 . . MAXNODE ? 

I D X = 1 . . M A X I D X ? 

RECLEN=1 . . MAXREC* 

DANUM-1. .MAX DA? 

FLAG=0 . .7? 

ADDR= INTEGER? 

POS I NT = 1 . .MAX I NT? 


<* PASLIB EXTERNAL PROCEDURES AND FUNCTIONS FOLLOW *> 

( * * ) 

<* TIMER PROCEDURES *> 

FUNCT I ON XT I ME : LONG I NT ? EXTERNAL ? 

PROCEDURE TSTART ( T : POS I NT ) ? EXTERNAL ? 

FUNCTION TREAD? LONG I NT? EXTERNAL! 

PROCEDURE TSTOP ! E X TERNAL ! 

<# TEXT OUTPUT PROCEDURES *> 

PROCEDURE MSG ( STR I NG : PACKED ARRAY C 1 . . ? 3 OF CHAR ) ! E X TERNAL ! 
PROCEDURE MSGLN ( STRING? PACKED ARRAY! 1 . . ?] OF CHAR ) ! EXTERNAL! 
PROCEDURE NXTLN? EXTERNAL! 

PROCEDURE MSG I < I : INTEGER) ! EXTERNAL! 

PROCEDURE MSGR < X : REAL ) ! EXTERNAL ! 

PROCEDURE MSGL < I .* LONG I NT ) ! E X TERNAL ! 

PROCEDURE ENDLN ( N : POS I NT ) ! EXTERNAL ! 

( ■»• * ) 

( * INTERACTIVE INPUT *> 


PROCEDURE QUERY! EXTERNAL! 

FUNCTION RDI: INTEGER! EXTERNAL! 

FUNCT I ON RDR : REAL ! EXTERNAL ! 

( # 

<3 DATA AREAS 

FUNCT I ON DAPTR < DA : DANUM ) : ADDR ! E X TERNAL \ 


<* FLOATING OPERATIONS *) 

FUNCTION ADD ( X , Y: REAL ) : REAL! EXTERNAL; 
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FUNCTION 

FUNCTION 

FUNCTION 

FUNCTION 

FUNCTION 

FUNCTION 

FUNCTION 

<* 

< * 

FUNCTION 

FUNCTION 

( * 

FUNCTION 

FUNCTION 

FUNCTION 

FUNCTION 


SUB ( X , Y : REAL > : REAL S EXTERNAL S 
MULT < X , Y : REAL ) : REAL S EXTERNAL 5 
DIVEK X, Y: REAL) : REAL* EXTERNAL; 
NEG< X: REAL ): REALS EXTERNALS 
ABS95 ( X : REAL ) : REAL ; EXTERNAL ; 
CMP < X , Y : REAL ) : REAL S EXTERNAL S 
SQRT95 < X : REAL ) : REAL; EXTERNAL; 


*) 

-*) 


— FLOATING POINT CONSTANTS 

MAX95: REAL; EXTERNAL; 

M I N95 : REAL S EXTERNAL; 

— FLOATING POINT CONVERSIONS *) 

I F I X ( X : RE AL ) : I NTEGER S E X TERN AL ; 

FLOAT I < I : INTEGER) : REAL; EXTERNAL; 

CV9512 < X : REAL ) : REAL; EXTERNAL ; 

VDP<N: POSINT; VAR A: ARRAY! 1 . . ?] OF REAL; 

VAR 8: ARRAY Cl. .?3 OF REAL >: REAL; EXTERNAL; 

FUNCTION S I NE<X: REAL): REAL; EXTERNAL; 

<* NBR COMMUNICATIONS *> 

PROCEDURE SEND2(N:N0DE; INDEX: IDX; LOC: ADDRS NWORDS: RECLEN) ; EXTERNAL; 
SEND2ALL ( INDEX : IDX ; LOC: ADDRS NWORDS: RECLEN ) ; EXTERNAL; 

RECV < N : NODE; LOC: ADDR ; NWORDS: RECLEN ) ; EXTERNALS 
RECV2<N:N0DE; INDEX: IDX; LOC : ADDR S NWORDS: RECLEN ); EXTERNAL; 
FLAG ROUTINES *) 


PROCEDURE 
PROCEDURE 
PROCEDURE 
<* 


PROCEDURE FLGEN ( F: FLAG ) ; EXTERNAL S 
PROCEDURE FLGSET ( F: FLAG ) ; EXTERNAL; 
PROCEDURE FLORES < F: FLAG ) S EXTERNAL; 
FUNCTION ANY (F: FLAG) : BOOLEAN; EXTERNALS 
PROCEDURE BAR (F: FLAG) ; EXTERNAL; 

(* PROCESSOR IDENTIFICATION 

FUNCTION PSELF: NODES EXTERNALS 
FUNCTION LSELF: NODES EXTERNALS 


-*) 


PROCEDURE P8CL**S 


CONST 

EDIM=3S 

TYPE 

TANDIM=1. .EDIMS 

TANMAT=ARRAYCTANDIM,TANDIM3 OF REAL; 
VECT0R2=ARRAY C T AND I M 3 OF REAL; 

DA4= ARRAY Cl. .91 OF INTEGER; 

DA5=ARRAY Cl. . 163 OF REALS 

VAR 
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S : T ANMAT ? EL : BOOLEAN ? 

F’B: REAL* NENTRY: INTEGER? DATA: ARRAYC 1 .. 43 OF REAL? 

NEW? 1 0 , MD I M , ND I M , NSEC , NPR 

CLK8.CLK9 
STRES, LOAD? LING, LOADS ? DEL 
PP, PR, MXP, MYP? MXRE, MYRE 
ZERO , ONE , TWO , FOUR ? TWLVE ? SIXTN? PI 
AR , ARND ? IX, I XND, I Y, I YND? ZX , ZXND, ZY, ZYND, RXND? RYND 
EBW, EBT, EDW, EAB, EAD, EDT , XRTB, XRTD, YRCB? YRCD 
IDX1 ? IDX2, I PR, NPROB, NCNTR, NS, NB, ND, NTB, NTD, CNTR, CNT2 
B, D, T, E, 3IGY, SIGRG, SIGRT, SIGRTM, SIGR, ER 

DET , RC , RT 

U I NT , V I NT , LENGTH , KBX, KBY, KTX , KT Y , SEGL 
<* Sprins constants #) CBX, CBY, CTX, CTY, Cl , C2 

DALL4 : @DA4 ? DALL5 
TOLi , T0L2 , T0L3 , T0L4 

<* SIGRTM MODIFIED TO a/c FOR c/c CALCULATIONS *> 
PROCEDURE PROPERTY*, 


ZBAR: REAL ? 
INTEGER? 
LONG I NT? 
VECTORS? 
REAL? 

REAL? 

REAL? 

REAL? 

INTEGER? 

REAL? 

REAL? 

REAL? 

REAL? 

®DA5? 

REAL? 


VAR 

B_2T , B_T , D_2T , D_T , BSQ , BCU , TSQ , TCU , DSQ , DCU : REAL ? 

FNB , FNTB , FND , FNTD : REAL ? 

BEGIN 

<* *) 
<* CALCULATION OF GEOMETRIC PROPERTIES *> 


<* CALCULATION OF AREA S< AREA ND *> 

AR : -ADD ( MULT ( TWO , MULT ( D , T ) ) , MULT ( TWO , MULT ( T , SUB ( B , MULT ( T , TWO )))))? 
ARND : =MULT ( FOUR , D I VD < AR , MULT ( B, D ) ) ) ? 

( *MSG( 'AREA OF CROSS-SECTION AR= ' > ? 

MSGR ( AR ) ? MSGLN ( ' I N**2 ' ) ? 

MSG ( 'NON-DIMENSIONAL AREA ARND-')? 

MSGR ( ARND ) ? ENDLN < 2 ) ? * ) • 

B..2T : —SUB ( B , MULT < TWO , T ) ) ? 

B._T : —SUB ( B , T ) ? 

D„2T : -SUB ( D , MULT ( TWO , T) ) ? 

D_T:=SUB(D,T>? 

TSQ : -MULT (T,T) ? 

TCU? -MULT ( TSQ, T ) ? 

BSQ: -MULT (B,B>? 

BCU: -MULT (B, BSQ) ? 

DSQ: -MULT ( D, D ) ? 

DCU: -MULT (D, DSQ) ? 
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( * CALCULATION OF IX ?< IXND 

I X : =D I VD ( MULT < B_2T , MULT ( T , MULT ( D_T , D_T ) ) ) , FOUR ) ; 

I X : =ADD ( IX,MULT(B_2T,DIVD(TCU,TWLVE> ) ); 

I X : = ADD (IX, MULT < T , D I VD ( DCU , TWLVE ) ) > ; 

I X : =MULT (IX, TWO ) ; 

( *MSG ( " MOMENT OF INERTIA I X= ' ) ; 

MSGR (IX); MSGLN ( ' IN**4'); *) 

I X NEi : =MULT ( B , DCU ) ; 

IXND: =MULT ( SIXTN, DI VD ( IX > IXND) ) ; 

(*MSG( "NON-DIMENSIONAL MOMENT OF INERTIA IXND* - > ; 
MSGR ( I XND ) ; ENDLN ( 2 ) ; * ) 

RXND: =SQRT95(DIVD( IXND, ARND) >; 

<* CALCULATION OF IY ?/ IYND 

I Y : =MULT ( D_2T , MULT ( T , MULT ( B_T , D I VD ( B_T , FOUR ) > ) ) ; 

I Y : =ADD ( IY,MULT(D_2T,DIVD(TCU, TWLVE) ) ); 

IY:=ADD( I Y , MULT ( BCU , D I VD ( T, TWLVE > ) ) ? 

IY: “MULT (TWO, IY); 

I YND : *MULT ( S I XTN , D I VD ( I Y , MULT ( BCU , D ) ) >; 

( *MSG( •'MOMENT OF INERTIA IY= '); 

MSGR ( I Y > ; MSGLN ( y I N**4 ' ) ; 

MSG ( 'NON-DIMENSIONAL MOMENT OF INERTIA I YND= •' ) ; 
MSGR ( IYND) ?ENDLN(2); *) 

RYND:=SQRT95(DIVD( IYND, ARND) ) ; 


( * PLASTIC MODULI ZX , ZXND, ZY , ZYND *) 

Z X : = ADD ( MULT ( B , MULT ( T , D I VD ( D_T , T WO ) ) ) , MULT ( T , MULT ( B_2T , 

D I VEi ( D_2T , FOUR ) )))' 

ZX:=MIJLT (TWO, ZX ) * 

(»MSG( "PLASTIC MODULUS ZX= ' ) ; 

MSGR ( Z X ) ; MSGLN ( ' IN**3 •' ) * *) 

Z X ND : =MULT ( F OUR , MULT ( TWO , D I VD ( Z X , MULT ( B , DSQ ) ) ) ) ; 

( *MSG ( •' NON-D I MENS I ONAL PL AST I C MODULUS Z XND= •' ) ; 

MSGR (ZXND); ENDLN (2);*) 

ZY: =ADD(MIJLT (D, MULT (T , DIVD(B_T, TWO) ) ) , 

MULT ( T , MULT ( B_2T , D I VD ( B_2T , FOUR ) ) ) ) ; 

ZY : =MIJLT ( TWO, ZY ) ; 

( *MSG ( - PL AST I C MODULUS Z Y= ' ) ; 

MSGR ( ZY ) ; MSGLN ( ' I N**3-' ) ; * ) 

ZYND : =MULT ( FOUR , MULT ( TWO , D I VD ( Z Y , MULT ( D , BSQ ) ) ) ) ; 

(*MSG( •'NON-DIMENSIONAL PLASTIC MODULUS ZYND* " ) ; 

MSGR ( Z YND ) 5 ENDLN ( 2 ) ; * ) 


(* 


ELEMENTAL PROPERTIES 


*) 
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<* ALL QUANTITIES ARE NON-DIMENSIONAL *) 

FNB ! =FLOAT I ( NB ) ? 

FNTB : =FLOAT I ( NTB ) ? 

FND:=FLOATI (ND> ? 

FNTD : -FLOAT I ( NTD ) ? 

EBW ! =D I VD ( TWO > FNB ) ? 

EBT : =MULT ( TWO, DIVD ( T i MULT (Dj FNTB) ) )? 

EDW : -MULT ( TWO > DIVB(T» MULT ( B , FNTD > > > ? 

EDT : =D I VD ( D_2T , FND ) ; EDT : =D I VD < MULT ( TWO , EDT ) , D ) 5 
<*M8GLN( 'MAJOR AXIS - HORIZONTAL' ) ? 

MSGLNt 'ELEMENTAL PROPERTIES - NOND I MENS I ONAL ' ) ? 

MSGLN < ' > • 

MSG ('WIDTH OF ELEMENT IN FLANGES EBW- ')? 

MSGR ( EBW ) * NXTLN ; 

MSG ('DEPTH OF ELEMENT IN FLANGES EBT= ' ) ? 

MSGR ( EBT ); NXTLN? 

MSG ('WIDTH OF ELEMENT IN WEB EDW- ')? 

MSGR ( EDW )? NXTLN? 

MSG ('DEPTH OF ELEMENT IN WEB EDT- ')? 

MSGR ( EDT > ? ENDLN ( 3 > ? * > 

IF CMP (RT, ZERO) =0 
THEN 

BEGIN 

XRTB: -ZERO? 

XRTD: =ZERO? 

YRCB? -ZERO? 

YRCD: =ZERO? 

END 

ELSE 

BEGIN 

XRTB: =ADD( B? D_2T ) ? 

XRTB : =MULT ( MULT ( RC , RT > , XRTB > ? 

XRTB : =D I VD ( XRTB , MULT ( MULT ( TWO , ADD ( RC , RT ) ) , ADD ( RC , RT ) ) ) ? 

XRTD: -MULT ( TWO » DIVD( XRTB» D) ) 5 
XRTB : -MULT ( TWO , D I VD ( XRTB , B ) ) ? 

YRCB : -MULT ( D I VD ( RC , RT ) , XRTB ) ? 

YRCD : -MULT ( D I VD ( RC , RT ) , XRTD ) ? 

S I GRTM : =D I VD ( MULT ( RT , T ) , MULT ( D , XRTD ) ) ? 

S I GRTM : -ADD ( RT , S I GRTM ) ? 

( * MSG ('SI GRTM- ' ) ? MSGR ( S I GRTM ) ? NXTLN ? 

MSG ('SI GRC- ' ) ? MSGR ( S I GRC ) ? ENDLN <3)5 * ) 

END? 

( *MSGLN ( ' RES I DUAL STRESS D I STR I BUT I ON D I MENS I ONS ' ) ? 

MSG ('IN FLANGES XRTB- ' ) ? MSGR ( XRTB > ? MSG ( ' YRCB- ' ) ? MSGR ( YRCB ) ? NXTI N ? 
MSG ('IN WEB XRTD- ')? MSGR ( XRTD )? MSG ( ' YRCD- ')? MSGR ( YRCD )? NXTLN? * ) 
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(* CALCULATE CONSTANTS *> 

SEGL : =FLOAT I ( NSEC- 1 > ; 

SEGL : =D I VD ( LENGTH , SEGL > ; 

C 1 : =MULT ( DSQ ? DIVD(E? SIGY) ) * 

C 1 : =D I VD ( C 1 ? MULT ( MULT ( SEGL ? SEGL ) ? FOUR ) ) ; (* Relates phi-ubar, vbar *) 

C2 : =MULT ( BSQ ? DIVD(E? SIGY) ) * 

C2 : =D I VD ( C2 , MULT ( MULT ( SEGL , SEGL ) , FOUR ) ) ; 

(* CBX? CBY? CTX? CTY — Constants of Restrints *) 

CBX : -MULT ( MULT ( TWO ? FOUR ) ? MULT ( I X , S I GY > > ; 

CBX : =MULT ( CBX ? SEGL > 5 CBY : =CBX ; 

CBX : =D I VD ( DSQ , CBX ) ; 

CBY : =D I VD ( BSQ ? CBY ) ? 

CTX : =MULT < CBX , KTX ) ? 

CBX : =MULT (CBX? KBX ) ; 

CT Y : =MULT ( CBY ? KTY ) ? 

CBY : =MULT ( CBY , KBY ) ? 

NXTLN ? MSG ( ' CBX ? CTX ? CBY ? CTY ? C 1 ? C2 are '); 

MSGR ( CBX ) ? MSGR ( CTX ) ; MSGR ( CBY ) ; MSGR ( CTY ) ? MSGR (Cl>! MSGR ( C2 ) ; NXTLN ; 

END*, (* PROCEDURE PROPERTIES *) 

( * # ) 

PROCEDURE ROUT I NE ( N ? M : I NTEGER ) *, 

LABEL 10; 

VAR 

COEF : ARRAY! 1 . . 143 OF REAL; I DSEC? ITER? TAN I DX: INTEGER; 

STORED , STOREL s ARRAY Cl. . NI OF REAL ; 

VCTRD ? VCTRL : ARRAY C 1 . . N ] OF RE AL ; 

UV s ARRAY C 1 . . M3 OF RE AL ; 

UO? VO: REAL; 

ITR: ARRAYC 1 . . NPR? 1 . . NPROB3 OF INTEGER; 

<* PROCEDURE TO CALCULATE RESIDUAL STRESSES *> 

PROCEDURE RESY ( VAR SIGR, EPSR: REAL; Y: REAL) ; 

VAR 

YD , FNB ? FNTB , XRD : REAL ; 

BEGIN 

XRD : = ADD ( XRTD ? D I VD ( T , D ) ) ; (* a/c FOR CORNER DEDUCTION #) 

YD:*ABS95(Y>; 

I F CMP ( YD ? SUB ( ONE , ADD (XRD, YRCD ) ) ) <=0 THEN 
BEGIN 

SIGR: =SIGRC; 

EPSR: =S I OR; 

END (* IF - THEN *> 

ELSE 
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( * 

end; 

PROCEDURE RESX ( VAR SIGR, EPSR: REAL; X: REAL) ; 

VAR 

XRB, XB,TT,T1:REAL; 

BEGIN 

XRB: “ADD( XRTB, DIVD(T » B) ) ; (* a/c FOR CORNER *> 

XB:=ABS95(X) ; 

I F CMP ( XB * SUB ( ONE , ADD ( XRB , YRCB ) ) ) OO THEN 
BEGIN 

SIGR: “SIGRC; 

EPSR: “SIGR; 

END <* IF - THEN *) 

ELSE 

BEGIN 

TT:=SUB(ONE,XB>; 

Tl:=ABS95< SIGRC); 

IF CMP (RT, ZERO) OO 
THEN 
BEGIN 

S I GR : =D I VD ( ADD (T1 » SI GRTM ) , ADD < XRB , YRCB ) ) ; 

S I GR : =SU8 ( S I GRTM , MULT < TT , S I GR ) ) ; 

END 

ELSE SIGR: = ZERO; 

EPSR: »SIGR; 

END; <* IF - ELSE *> 

( # MSG ( •' RES X - ) ; MSGR ( S I GR ) ; MSGR ( EPSR ) ; N X TLN ; * ) 

END; <* PROCEDURE END *) 

PROCEDURE CHECK ( VAR EL : BOOLEAN ) ; FORWARD ; 

PROCEDURE ESTIFF (VAR S: TANMAT ; VAR EL: BOOLEAN; DEL: VECTORS; 


BEGIN 

FNB : =SUB ( ONE , YD ) ; 

FNTB : = ABS95 ( S I GRC ) ; 

IF CMP (RT, ZERO) OO 
THEN 
BEGIN 

S I GR : =D I VD ( ADD ( FNTB , S I GRTM ) , ADD ( X RD , YRCD ) ) ; 
S I GR : =SUB ( S I GRTM , MULT ( S I GR , FNB ) ) ; 

END 

ELSE SIGR: “ZERO; 

EPSR: “SIGR; 

END; (* IF - ELSE *) 

MSG ( •" RES' Y •" ) ; MSGR (SIGR); MSGR < EPSR ) ? NX TLN ? * ) 

(* PROCEDURE END *) 
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ARND, IXND, IYND: REAL' VAR PP, PR, MXP , MYP, MXRE, MYRE: REAL ); FORWARD? 
PROCEDURE SOLVE ( X : TANMAT ; SUMS : VECTORS 5 VAR SOLNS : VECTORS ; 

VAR DET : REAL ) ; FORWARD ; 

PROCEDURE TANGNT(LD, LOADS: VECTORS; VAR DEL: VECTORS; VAR S: TANMAT; 

VAR CNT2: INTEGER; UO, VO: REAL); 

LABEL 10; 

VAR 

LD 1 , SOLNS , DEL 1 : VECTORS ; 

CNT : INTEGER; 

CLK 1 , CLK2 : LONG I NT ; 

BEGIN 

cntss-o; 

CHECK ( EL ) ; 

IF (EL) THEN 
BEGIN 

SOLVE ( S , LOADS , SOLNS , DET ) ; 

FOR I : = 1 TO EDIM DO DEL C I 3 : =SOLNS C I 3 ; 

END (* — IF - THEN — *) 

ELSE 


BEGIN 

FOR I : =1 TO EDIM DO DELI C I 3 : =DELC I 3 *, 

REPEAT <* — UNTIL LDC1,2,&33 < TOL — *> 
CNT : =0 ; 


CNT2:=CNT2+l; 


FOR I : =1 TO EDIM DO LD C I 3 : =SIJB ( LOADB C I 3 , LD C I 3 ) ; 


FOR I : =1 TO EDIM DO 

t p CMP ( ABS95 ( LD Cl 3 ) , TOL 
IF CNT<EDIM "THEN 
BEGIN 


— /“\ t I i 1“ k i r-iiiT . — .-.nr . 4 . 
^ V I I m.M I ■ H I * X 7 


SOLVE ( S , LD , SOLNS , DET ) ; 

IF CMP ( DET » T0L2 ) <=0 THEN GOTO 10; 

FOR I : =1 TO EDIM DO DELI C I 3 : =ADD ( DELI C I 3 , S0LNSCI3 ) ; 

CALCULATE STRESSES FI *> 

ESTIFF(S,EL,DEL1, ARND, IXND, I YND, PP, PR, MXP, MYP, MXRE, MYRE) ; 
FOR I i=l TO EDIM DO LD1 C I 3 : =VDP ( EDIM, SC I 3 , DELI ) ; 

LD1 C 1 3 : =ADD ( LD 1 C 1 3 , ADD ( PP, PR ) ) ; 

LD1 C23 : =ADD(LD1 C23 , ADD (MXP, MXRE) ) ; 

LD 1 C 3 3 : =SUB ( LD 1 C 3 3 , ADD ( M YP , MYRE ) ) ; 

FOR I : =1 TO EDIM DO LD C I 3 : =LD 1 C I 3 ; 

END; (* ELSE CNVRGOTRUE *> 

IF CNT=EDIM THEN 

FOR I : =1 TO EDIM DO DELL I 3 : =DEL1 C I 3 ; 

IF CNT2=15 THEN MSGLN ( '' NO CONVERGENCE IN 15 ITERATIONS ' ) 
UNTIL (CNT=EDIM) OR (CNTS>=15) 

END; (* ELSE *> 
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<* MSGLN ( - CONVERGED VALUES OF LOAD AND DISPLACEMENT ' ) ; 

MSGLN ( ' ) ? 

MSG < •" •" ) ; MSG ( - LOAD VECTOR ' ) ; 

MSG ( •' - ) ; MSGLN ( " DISPLACEMENT " > ; 

FOR I s =1 TO EDIM DO 
BEGIN 

NX TLN ? MSG Id); MSG (' ')% MSGR ( LOADB C I 3 > 5 MSG ( ' x ) ; MSGR ( DEL EH)? 

END; *) 

10: END; <* PROCEDURE TANGNT *) 


<# PROCEDURE ESTIFF 

PROCEDURE ESTIFF; (* — (VAR S : T ANM AT ; VAR EL: BOOLEAN; 

ARND, IXND, I YND: REAL; VAR PP,PR,MXP,MYP,MXRE,MYRE:REAL>- 

VAR 

CNT 1 , CNT2 , MB , MD : I NTEGER ; 

SUM , ASUM , BSUM , CSUM , DSUM , ESUM , T X , T Y , X , Y , TMP : REAL ; 

DEL2, DELS, B: REAL; 

PROCEDURE STIFF; 

NO OPTIMIZE *) 

BEGIN 

SC 1 , 1 ] : =DI VEKSUM, ARND) ; 

SC 1,21: =DIVD( ASUM , ARND ) ; 

PP:=DIVD(PP, ARND); 

PR : =D I VD ( PR , ARND ) ; 

MXP:=DIVD(MXP, IXND); 

MYP : =D I VD ( MYP , I YND ) ; 

MXRE: =DIVD(MXRE, IXND); 

MYRE : =D I VD ( MYRE , I YND ) ; 


S C 1 , 3 3 : =NEG ( D I VD ( BSUM , ARND ) ) ; 

SC 3, 13: =NEG ( D I VD < BSUM , I YND ) ) ; 

S C 3 , 2 3 : =NEG ( D I VD ( CSUM , I YND ) ); 

SC3, 33 : =DI VD< DSUM, I YND ) ; 

S C 2 , 3 3 : =NEG ( Ci I VD ( CSUM , IXND) ); 

S C 2 , 2 3 : =D I VD ( ESUM , IXND); 

SC2, 1 3 : =Ei I VD < ASUM , IXND) ; 

end; 

PROCEDURE ST I FCAL ( MB , NTB , I DX : I NTEGER ; 

EBT , EBW , EAB , DEL2 , DEL3 : REAL ) ; 


*> 

*> 


LABEL 20,30; 


BEGIN 

<* IDX=0 USES RESX IDX=1 USES RESY *> 

CNT 1 : = 1 ; 

30: FOR I : =1 TO NTB DO 

BEGIN 
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Y:=FLOATI (2*1-1 ) ; 

Y : =SUB ( ONE , MIJLT < Y , D I VD ( EBT , TWO > > > ; 

IF CNT 1 >=2 THEN Y : =NEG ( Y ) ? 

TY:=MULT(Y,DEL2); 

CNT 2 s =1 ; 

20: FOR J : - 1 TO MB DO 

BEGIN 

X : =FLO AT I < 2* J- 1 ) * 

X:=MULT(X,DIVD<EBW,TWO> ) ; 

IF CNT2>=2 THEN X:=NEG(X); 

IF IDX=0 THEN RESX (SIGR, ER, X ) 

ELSE RES Y ( S I GR » ER , X ) 5 
TX:=MULT(DEL3, X) ; 

IF IDX=0 THEN TX : =ADD ( SUB ( DELC 1 ] , TX ) , ADD ( TY , ER ) ) 
ELSE TX : =ADD ( SUB ( DELC 1 I , TY > , ADD < TX , ER > > ; 
I F CMP ( ABS95 < TX > , ONE ) >=0 THEN 
BEGIN 

IF CMP (TX, ZERO )<0 
THEN 
BEGIN 

PP : =ABD ( PP , NEG ( E AB ) ) ; 

IF IDX=0 THEN 
BEGIN 

MXP : =ADD (MXP, MULT (Y, NEG (EAB) ) ) ; 

M YP : = ADD < MYP , MULT ( X , NEG ( E AB ) > ) ? 

END 

ELSE BEGIN 

MXP: =ADD(MXP, MULT (X, NEG (EAB) ) ) ; 

MYP: =ADD ( MYPi MULT ( Y » NEG ( EAB > > ) ? 

END; 

END 

ELSE 

BEGIN 

PP : =ADD ( PP > EAB ) ; 

IF IDX=0 THEN 
BEGIN 

MXP:=ADD(MXF*,MULT(Y,EAB) ) ; 

MYP: =ADD(MYPi MULT ( X , EAB ) ) ; 

END 

ELSE BEGIN 

MXP : =ADD (MXP, MULT (X, EAB) ) ; 

MYP: =ADD( MYP, MULT (Y, EAB) ) ; 

END; 

END; <* — IF ELSE — *) 

END <* IF TX>1.0 *) 
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ELSE 

BEGIN 

SUM: =ADD(SUM, EAB) ? 

PR : = ADD ( PR , MULT (SIGR, EAB) ) ? 

IF IDX=0 THEN 
BEGIN 

ASUM : =ADD ( ASUM , MULT ( Y , EAB ) ) ? 

BSUM : = ADD ( BSUM , MULT ( X , EAB ) ) ? 

CSUM : -ADD ( CSUM , MULT ( MULT ( X , Y ) , EAB ) ) ? 

DSUM: =ADD ( DSUM , MULT ( MULT ( X , X ) , EAB > > ? 

ESUM : -ADD ( ESUM , MULT ( MULT ( Y , Y ) , E AB ) ) ? 

MXRE : -ADD ( MXRE > MULT ( S I GR , MULT ( Y , EAB ) ) > ? 

MYRE : -ADD ( MYRE , MULT ( S I GR , MULT ( X, EAB) ) ) * 

END 

ELSE BEGIN 

ASUM: -ADD ( ASUM , MULT ( X , EAB ) ) ? 

BSUM: » ADD ( BSUM, MULT ( Y, EAB ) ) ? 

CSUM: =ADD ( CSUM , MULT ( MULT ( X, Y) , EAB ) ) ? 
DSUM: =ADD ( DSUM , MULT ( MULT (Y, Y) , EAB ) ) ; 
ESUM: = ADD ( ESUM , MULT ( MULT ( X , X ) , EAB ) ) ; 
MXRE: =ADD ( MXRE , MULT ( S I GR , MULT ( X , EAB ) > ) ; 
M YRE := ADD ( M YRE, MULT ( SI GR, MULT (Y, EAB) ) ) ? 
END? 

END? <* — ELSE — *> 

END? <* FOR J *> 

CNT2: =CNT2+1 ? 

IF CNT2<=2 THEN GOTO 20? 

END? (* FOR I *) 

CNT 1 : =CNT 1 + 1 ? 

IF CNT 1 <=2 THEN GOTO 30? 

END? 

BEGIN 

FOR I : =1 TO EDIM DO 
BEGIN 

FOR J: =1 TO EDIM DO SCI , JD : -ZERO 
END? 

CHECK (EL)? (* CHECK YIELDING OF CROSS-SECTION *) 

IF NOT (EL) THEN 
BEGIN 

MB: =NB DIV 2? 

MD:=ND DIV 2? 

PR: -ZERO? 

PP: =PR; 

MXP: =PP? 

MYP: =PP? 
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MX RE: =PP; 

MYRE: -PP; 

<* note: P -VE, MX +VE, and MY -VE 

phiY BE P03. VALUE & phiX BE POS. VALUE *> 

EAB: -MULT ( EBW, EBT) ; 

EAD s -MULT ( EBW » EOT > * 


!r« 

UM : = 

ZERO 

A 

SUM: 

-SUM 

B 

SUM: 

-SUM 

c 

SUM: 

-SUM 

D 

SUM: 

-SUM 

E 

SUM: 

-SUM 


ST I FCAL ( MB , NTB , 0 , EBT , EBW , E AB , DEL C 2 3 , DEL C 3 3)5 
ST I FCAL ( MD , NTD , 1 , EDW , EDT , EAD , DEL LSI, DEL C 2 ] ) ? 
STIFF! 

END 

ELSE (* IF ELASTIC *> 

BEGIN 

FOR I : =1 TO EDIM DO 
BEGIN 

FOR J:=l TO EDIM DO 
BEGIN 

IF IOJ THEN SCI, J]: -ZERO 
ELSE S C I , J3 : —ONE ‘ 

END; 

END; 

PR: -ZERO; 

PP: =F’R; 

MXP: =PP; 

MYP : =PP ; 

MXRE: =PP; 

MYRE : -PP ; 

END; <* — ELSE -- #> 

END; (* — PROCEDURE ESTIFF -- *) 


<* PROCEDURE TO CHECK YIELDING OF CROSS-SECTION *) 

<* *> 


PROCEDURE CHECK; <* —(VAR EL: BOOLEAN) — *> 
VAR 

X » Y? XRB, XRD: REAL; 
is INTEGER; 
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BEGIN 

XRBs =AOD< XRTB> DIVD(T> B) )? 
XRD:=ADD<XRTD,DIVD<T,D> )? 

X : “ADD < XRB » YRCB ) * 

X s “SUB ( ONE iX)i 
X : =D I VD < X j EBW ) ’ 

I ! =IFIX ( X ) * 

I : =2# ( 1-1 ) + l ; 

X s =FLOAT I < I ) 5 
X : “MULT < X , D I VD < EBW » TWO ) ) ? 

IF CMP < DEL C 3 3, ZERO) <0 THEN X:=NEG<X)5 
RESX (SIGRj ERj X ) * 

X i =MULT ( X 7 DEL u 3 j > » 

Y : “SUB ( ONE , D I VD ( EBT , TWO ) ) ; 

IF CMP < DEL C 2 3, ZERO) >0 THEN Y:=NEG<Y); 

Y : “MULT ( Y ? DEL C 2 3 ) * 

X s = ADD ( SUB ( DEL C 1 3 > X ) ? ADD < Y » ER > ) ? 

IF CMP ( ABS95 ( X ) » ONE ) >=0 
THEN EL: =FALSE 
ELSE 
BEGIN 

Y:=ADD(XRD,YRCD); 

Y : =SUB ( ONE » Y ) * 

Y s =D I VD ( Y » EDT ) * 

I : =IFI X ( Y ) ’ 

I : =2# i X — i ) + i; 

Y: =FLOATI ( I ) i 
Y : =MULT ( Y , D I VD ( EDT , TWO > ) ; 

IF CMP ( DEL C 2 3 » Z ERO ) >0 THEN Y : =NEG < Y ) * 
RESY(SIGR,ER, Y) ? 

Y : =MULT ( Y ■> DEL C 2 1 ) * 

X : =SUB ( ONE , D I VD ( EDW , TWO ) ) ; 

IF CMP < DEL C 3 I, ZERO XO THEN X:=NEG(X); 

X : =MULT ( X ? DEL C 3 3 ) ; 

Y : =ADD ( SUB ( DEL C 1 3 > X ) * ADD < Y * ER ) ) ; 

IF CMP ( ABS95 ( Y ) > ONE ) >=0 

THEN EL: “FALSE 

ELSE 

BEGIN 

X : =SUB ( ONE , D I VD ( EBW , TWO ) ) ; 

I F CMP ( DEL C 3 1 , Z ERO ) >0 THEN X : =NEG ( X ) ; 
Y : =SUB < ONE , D I VD ( EDT , TWO ) ) ; 

IF CMP ( DEL C 2], ZERO) <0 THEN Y : =NEG ( Y ) * 
X : “MULT ( X i DEL C 3 ] ) * 

Y : “MULT ( Y •» DEL C 2 3 ) ? 


C - oL 
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X : =ADD ( SUB ( DEL C 1 3 , X ) , ADD < Y> SI GRT ) ) ; 

IF CMP ( ABS95 ( X ) > ONE ) >=0 
THEN EL:=FALSE 
ELSE EL:=TRUE; 

END? <* — ELSE — *> 

END; <* — ELSE — *> 

IF NOT EL THEN FLGSET ( 5 ) ; 

END; (* PROCEDURE *) 

<* MATRIX INVERSION SUITABLE FOR FEM *) 

( w#*#*********###**********-**#*#--*******-###*#*-#-#-*-*##*###*###-* ) 

PROCEDURE MATINV ( VAR AsTANMAT; N: TANDIM) ; 

<* MATRIX INVERSION PROCEDURE *> 

<* METHOD - JOHNSON ( 1 ) METHOD *> 

<* A IS N X N SYMMETRIC MATRIX*) 

VAR 

NSWAP: INTEGER; IPV: ARR A YC TANDIM, 1 . . 33 OF INTEGER; 

L» I * I COL 7 JCOL , I ROW 7 JROW : T AND I M ; 

T 7 P I V , SWAP , AMAX : REAL ; 

BEGIN 

DET : =ONE; 

FOR J: =1 TO N DO IPVC J, 33 : =0; 
i:«l; 

REPEAT (* UNTIL DET =0 OR I>N *) 

AMAX:=ZERO; 

FOR j:=l TO N DO 
BEGIN 

IF I P V C J , 3 3 0 1 THEN 
BEGIN 

FOR K : = 1 TO N DO 
BEGIN 

IF IPVCK7 3301 THEN 

IF CMP ( AMAX 7 ABS95 ( AL J, K3 ) ) <0 THEN 
BEGIN 

IROW:=J; (* SAVE ROW ?< COLUMN *) 

icol:=k; 

AMAX : =ABS95 ( AC J , K3 ) ; 

END; (* IF AMAX *> 

END; <* FOR K *) 

END; (* IF IPVCJ733 *) 

END; (* FOR J *) 

(* IF PIVOT ELEMENT IS EQ 0 - DET = 0 *) 
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IF CMP ( AMAX , T0L2 ) <=0 THEN DETs “ZERO 
ELSE 
BEGIN 

I P V C I COL , 3 ] : = I P V C I COL , 3 ] + 1 ; 

I PV Cl, 1 ] s = I ROW ' 

I P V C 1 , 2 ] ; = I COL * 

(* INTERCHANGE ROWS TO PUT PIVOT ELEMENT ON 
IF I ROMO I COL THEN 

FOR J: =1 TO N DO 
BEGIN 

SWAP s = A n I ROW iJl* 

AC I ROW, J] s =AC ICOL, J] ; 

A C I COL » J ] s “SWAP ; 

END? (* FOR J 1 TO N *) 

<* DIVIDE PIVOT ROW BY PIVOT ELEMENT *> 
PIV:=ACICOL, ICOL]; 

DET ! =MULT ( DET » P I V ) ; 

AC ICOL, ICOL] s “ONE* 

FOR K:=l TO N DO 

A C I COL , K ] : =D I VD ( A C I COL » K ] , P I V ) ; 

<* REDUCE THE NON-PIVOT ROW BY SUBTRACTION *> 
FOR LI : =1 TO N DO 

IF LIOICOL THEN (* AVOID DIAGONAL *) 
BEGIN 

T:=ACL1, ICOLl; 

A C L 1 , I COL ] s = ZERO ; 

FOR K:=l TO N DO 

A C L 1 » K ] : =SUB< ACL1 , K] , MULT ( AC ICOL, KD , T ) ) 
END; <* IF LI O ICOL *) 

I ; =1+1 ; (* GO ON TO NEXT *) 

END; <* ELSE IF AMAXCTOL *) 

UNTIL ( I >N ) OR CMP (DET, ZERO )=0; 

<* INTERCHANGE THE COLUMNS AND MODIFY DET *> 
NSWAPs =0; (* THIS IS SIGN FLAG *) 

IF CMP (DET, ZERO) OO THEN BEGIN (# DO ONLY 
FOR j: =1 TO N DO 
BEGIN 
Ls-N-J+l; 

IF I F'V C L , 1 ] <> I P V C L , 2 ] THEN 
BEGIN 

JROW ; = I P V CL, 1 ] ; 

•JCOL : = I P V C L , 2 ] ; 

NSWAP: =NSWAF'+1 ; (* COUNT SWAPS *) 

FOR K:=l TO N DO 
BEGIN 


DIAGONAL *) 


IF DET. NE. 0.0 *> 
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SWAP: =ACK, JR0W3S 
A C K , JROW 1 ! = A C K , JCOL 1 ; 

ACK, .JCOL3 : =SWAP; 

END? <* FOR K *) 

END; (* IF IPV *> 

END; <* FOR J *) 

IF ODD ( NSWAP > THEN SWAP: =NEG( ONE) 

ELSE SWAPS »ONE; 

END; (# IF DETOO.O *> 

DET : -MULT ( DET , SWAP ) ; 

END; (* MAT I N V PROCEDURE *) 

PROCEDURE SOLVE; 

VAR SIZE: INTEGER; 

BEGIN (* MAIN PROGRAM *) 

SIZE: =EDIM; 

MATINV ( X » SIZE) * <* CALL MATINV PROCEDURE *) 

IF CMP ( DET , T0L2 ) <1=0 
THEN MSGLN< 'MATRIX IS SINGULAR') 

ELSE 

BEGIN 

FOR J : = 1 TO SIZE DO 
BEGIN 

SOLNS C J I : = Z ERO ; (* INITIALIZE *) 

SOLNS C J 3 s =VDP ( SI ZE , X C J] , SUMS ) ; 

END; (* FOR J *) 

END; <* ELSE i.e. DET .NE. 0.0 *) 

END; 

BEGIN ( * Routine's main body starts here *) 

ITER: =0; FOR I:=l TO NPR DO FOR J:=l TO NF’ROB DO ITRC I , . I] : =D; 

U I NT : =D I VD < MULT ( U I NT , LENGTH > , D I VD ( B , TWO ) > ; 

VINT: =0 1 VD ( MULT (VINT, LENGTH ) , D I VD ( D , TWO ) ) ; 

REF’EAT (# Until failure of column *) 

I PR: =0; IDX 1 : =0; IDX2:=0; 

TANIDX: =0; 

CLKS: =TREAD; 

BAR ( 2 ) ; 

CLK9 : =CLK9+TREAD-CLK8 ; 

IF ANY (6) THEN GOTO 10; 

CLK8:=TREAD; 

RECV2(NPR+1, 1 , LOCATION <PB) , 11 ) ; (-B-Recieve PB to DAT A [ 4 3 ; DAT A=M coef*) 
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RECV2 ( NPR+ 1 1 2 , LOCAT I ON ( UV ) , 2*MD I M ) ; <*Recieve UBARs & VBARs *) 
CLK9 : = C: L K 9+TREAD-CL K 8 ; 

FOR I : =1 TO MDIM DO 

I F < ODD ( I ) ) THEN UV Cl]: =MULT ( D I VD ( PB , MULT ( RYND , RYND ) ) , UV C I ] ) 

ELSE UVC I 3 : =MULT < D I VD < PB, MULT ( RXND, RXND > ) , UVC I 3 ) ,* 
DATAC 1 3 : =MULT ( CBX , DATAC 13); DATAC 23 : =MULT ( CTX , DATAC 2 3 ) ; 

DATA C 3 3 : =MULT ( CBY, DATA C 3 3 ) ; DATA C 4 3 : =MULT ( CTY , DATA C 4 3 ) ; 

REPEAT (* Until NF'ROB are completed #) 

I PR: =IPR+1 * 

I TER : = ITR C LSELF , I PR 3 ; 

IF ( NENTRY=0 ) (* NENTRY=0 for start of GLOBAL PROBLEM *) 

THEN FOR I : =1 TO EDIM DO 

BEGIN DEL C I 3 : =ZERO ; LOADC I 3 : =ZERO END 
ELSE FOR I « =1 TO EDIM DO 

BEG I N DEL C I 3 : =STOREDC I + 1 DX 1 3 ; LOAD C I 3 : *STOREL C I + I DX 1 3 END ; 
EST I FF < S , EL , DEL , ARND , IXND, I YND, PP, PR, MXP, MYP, MXRE, MYRE ) ; 

IDSEC: =NPROB*( LSELF- 1 >+IPR; 

IF I DSEC= 1 THEN BEGIN UVC 1 3 : =ZERO; UVC23:=ZERO; END; 

IF I DSEC=NSEC THEN BEGIN UVCMDIM3 : =ZERO; UVCMDIM-1 3 : =ZERO; END; 
ZBAR! =DIVD<FLOATI ( IDSEC-1 ) , FLOAT I ( NSEC- 1 ) > ; 

UO : =MULT ( S I NE ( MULT (PI, ZBAR ) ) , U I NT ) ; 

VO: =MULT (SINE (MULT (PI , ZBAR) ) ,VINT>; 

LOADBC 1 3 : =NEG ( PB ) ; 

L0ADBC2 3 : =MULT ( SUB < ZBAR , ONE ) , DATAC 13); 

LOADS C 2 3 : =SUB ( ADD ( LOADBC23 , UV C IDX2+2 3 ) , MULT < ZBAR , DATA C 2 3 ) ) ; 
LOADBC 23 : =ADD ( MULT ( DI VD<PB, MULT < RXND, RXND ) ) , VO) , LOADBC 2 3 ) ; 

LOADS C 3 ] : =MULT ( SUB ( ONE , ZBAR ), DATA C33 ) ; 

LOADB C 3 3 : =ADD ( SUB ( LOADB C 3 3 , UV C I D X2+ 13), MULT ( ZBAR , DATA C 4 3 ) > ? 
LOADB C 3 3 : =SUB ( LOADB C 3 3 , MULT ( D I VD ( PB , MULT ( RYND , RYND ) ) , UO ) ) ; 

T ANGNT < LOAD , LOADB , DEL , S , CNT2 , UO , VO ) ; 

TANIDX : =0; 

ITRCLSELF , I PR 3 : = I TR C LSELF , IPR3 + 1; ITER: =ITRCLSELF, IPR3; 

IF 1003 AND ANY (5) THEN 
BEGIN NXTLN; FOR I : = 1 TO EDIM DO 

BEGIN FOR J: =1 TO EDIM DO MSGR(SC I , J3 ) ; NXTLN; END; END; 
IF 1004 AND ANY (5) THEN 

BEG I N MSG < x LOAD= ' ) ; MSGR ( PB ) ! MSG ( ' DET= y ) ; MSGR ( DET ) ; NXTLN ; END ; 

IF CMP ( DET , T0L2 ) OO OR CNT2=15 THEN TANIDX :=l; 

FOR I : =1 TO 14 DO COEFC 1 3 : =ZERO; 

IF TANIDX <> 1 THEN 
BEGIN 

COEF C S 3 : =SUB ( S C 2 , 3 3 , MULT ( S C 2 , 1 3 , DI VD ( SC 1 , 33 , SC 1 , 13))); 

COEF C 8 3 : =MULT ( C2 , COEF C 3 3 ) ; 

COEF C 9 3 : =MULT ( NEG < TWO ) , COEF C 8 3 ) ; 

COEFC 10 3 : =SUB ( MULT ( SC2, 1 3 , DI VD( SC 1 , 23 , SC 1 , 1 3 ) ) , SC2, 23 ) ; 
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COEF r. 1 0 3 : =MULT ( C 1 > COEF CIO])' 

COEFC m:=SUB (MULT (NEG( TWO ) , COEFC 103 ) , DIVD ( F’B, MULT ( RXND, RXND ) ) > ’ 
COEF C 12] : =MULT < CBX , SUB ( ZBAR , ONE ) ) 5 
COEF CIS]: =MULT ( CTX , NEG ( ZBAR ) ) ? 

COEF C 1 4 ] : =ADD ( PR , ADD ( PP , PB ) ) ; 

COEF C 1 4 ] : =ADD ( MULT ( PB , D I VD ( VO , MULT ( RXND , RXND ) ) ) , 

D I VD < MULT ( COEF C 1 4 ] , S C 2 , 1 ] ) , SC 1 , 1 ] ) ) ? 
COEF C 1 4 ] : =SUB < COEF C 1 4 ] , ADD ( M XRE , M XP ) ) ? 

FOR I : =8 TO 14 DO COEFC I ] : =NEG ( COEFC I ] ) ; <* Consistent with tan*) 

C0EFC1]:=SUB(SC3,3],MULT(SC3, 1 ] , DIVD(SC 1 , 3] , SC 1 , 1]) ) >; 

COEF Cl]: =MULT < C2 , COEF Cl])? 

COEF C 2 ] : =ADD ( MULT < NEG ( TWO ) , COEF C 1 ] ) , D I VD ( PB , MULT < RYND , RYND ) ) ) ; 

COEF C 3 ] : =SUB < MULT < S C 3 » 1 ] , DI VD ( SC 1 , 2] , SC 1 , 1 ] ) > , SC3, 21 ) i 

COEF C 3 ] : =MULT (Cl, COEF C 3 ] > 5 

COEF C 4 ] : =MULT ( NEG ( TWO ) , COEF C 3 ] ) ; 

COEF C 5 ] : =MULT ( CBY , SUB ( ONE , ZBAR ) ) ; 

COEF C 6 ] : =MULT ( CT Y , ZBAR ) ! 

COEF C 7 ] : =MULT ( D I VD ( S C 3 , 1 ] , S C 1 , 1 ] ) , ADD < ADD ( PR , PP ) , F'B ) ) ", 

COEF C 7 ] : = ADD ( ADD ( COEF C 7 ] , MYRE ) , SUB ( MYP, MULT ( PB , 

D I VD ( UO , MULT ( RYND , RYND ) ) ) ) ) ? 

FOR is-l TO EDIM DO 

BEGIN VCTRDC IDX 1 + 1 ] s =DELC I ] ? VCTRLC IDX 1 + 1 ] : =LOADBC I ] ; END; 

END; 

CLK8 : =TRE AD ; SEND2 ( NPR+ 1 , I PR , LOCAT I ON ( COEF ) , 3 1 ) ; 

CLK9 : -CLK9+TRE AD-CLK8 ; 

IDXl:»IDXl+EDIM; 

IDX2: =IDX2+(EDIM-1 ) ; 

UNTIL (IPR=NPROB) OR TANIDX=1; 

CLKS: =TREAD; BAR(7>; RECV2( NPR+1 , 3, LOCATION ( NEW ), 1 ) ; 

CLK9 : =CLK9+TRE AD-CLKS ; 

IF NEW=0 THEN 
FOR I : =1 TO N DO 

BEG I N STORED Cl]: =VCTRD Cl]; STOREL Cl]: =VCTRL C I ] END ; 

UNTIL ANY ( 6 ) ; <* Failure of column *) 

10:NXTLN; 

END; (* Routine Ends Here #) 

(* MAIN PROGRAM BEGINS *> 

BEGIN 

FLGEN ( 2 ) ; FLGEN ( 5 ) ; FLGEN ( 7 ) 5 FLGEN ( 6 ) ; 

FLGRES ( 2 > ; FLGRES ( 5 ) ; FLORES ( 7 ) ; FLGRES < 6 ) ; BAR ( SYSFLAG ) ; 

ZERO : =CV95 12(0.0); 

ONE: =CV9512( 1.0) ; 

TWO : =C V95 12(2.0) ; 

FOUR* -CV95 12 (4.0) ; 

T WL VE : =CV95 12(12.0); 
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SIXTN: =CV95 12(16.0? 

P I : =C V95 12(3.141 592654 ) ? 

<* GET DATA AREAS *> 

DALL4? : ADDRs =DAPTR( 4 > ? 

NF'ROB : =DALL4© CIO NCNTR : =DALL4@ C 2 3 ? 

NB: =DALL4©C33 ? ND? =DALL4©C43 ? 

NTB : =D ALL4© C 5 3 ? NTD : =D ALL4© C 6 3 ? 

NSEC: =DALL4@C73 ? NPR: =DALL4@C83 ? 10: =DALL4@C93 ? ( #NSEC-No. of Sections 

NPR -No. of Processors used less one #> 

DALL5: : ADDR : =DAPTR ( 5 ) ? 

B : =CV95 1 2 ( D ALLS® C 1 3 ) ? Da =CV95 1 2 ( BALLS© E 2 3 ) ? T : =CV95 1 2 ( D ALLS© C 3 3 ) ? 

RT : =C V95 1 2 ( BALLS© C 4 3 > ? RC : =CV95 1 2 ( BALLS© C 5 3 ) ? 

T0L2 : =CV95 1 2 ( BALLS© C 6 3 ) ? 

LENGTH: =CV9512(DALL5@C73 )? 

KB X : =CV95 1 2 ( BALLS© C 8 3 ) ? KB Y : =C V95 1 2 ( BALLS© C 9 3 ) ? 

KT X : =CV95 1 2 ( BALLS© C 1 0 3 ) ? KTY *. =C V95 1 2 ( BALLS© C 1 1 3 ) ? 

U I NT : =CV95 1 2 ( BALLS© C 1 2 3 ) ? VINT: =CV95 1 2 ( BALLS© C 1 3 3 ) ? 

E : =CV95 1 2 ( BALLS© C 1 4 3 ) ? S I QY a =CV95 1 2 ( BALLS© C 1 5 3 > ? 

TOL 1 : =C V95 1 2 ( BALLS© C 1 6 3 ) ? 

S I GRC : =NEG ( RC > ? 

S I GRT : =RT ? NXTLN? 

IF CMP (RT, ZERO) =0 THEN MSGLN C MEMBER WITHOUT RESIDUAL STRESSES'' ) 

ELSE MSGLN ('MEMBER WITH RESIDUAL STRESSES')? 

(*IF CMP (KBX, ZERO )=0 AND CMP (KTX, ZERO) =0 THEN 

MSGLN ('NO PARTIAL RESTRAINT AGAINST X AXIS)? 

IF CMP (CBY, ZERO )=0 AND CMP (CTY, ZERO) =0 THEN 

MSGLN ( x NO PARTIAL RESTRAINT AGAINST Y AXIS")? *) 
TSTART ( 1 ) ? CLK9: =0? 

PROPERTY? 

BET : =ONE? 

NPROB: =NSEC DIV NPR? 

NDIM:=EDIM*NPROB? 

MDIM: =(EDIM— 1 )*NPROB? (* NDIM, MDIM - Variable lengths of STORE, UV *) 
ROUT I NE ( ND I M , MD I M ) ? 

TSTOP? 

ENDLN ( 2 > ? MSG ( ' PROCESSING TIME =")? MSGL (TREAD) ? NXTLN? 

NXTLN? MSG( 'EXECUTION TIME =•' > ? MSGL ( TREAD-CLK9 ) ? NXTLN? 

END? (* End of the main program ■*) 

BEGIN <**N0 OBJECT * ) 

END . 
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APPENDIX D 


PROOF-COPY OF THE CONCURRENT PROCESSING PAPER FROM 
ENGINNERING WITH COMPUTERS JOURNAL 


A proof-copy of the paper mentioned in Section 3.3 is given in this 
appendix. 
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APPENDIX E 


FRAME SUBSTRUCTURING 


E. 1 Introduction 

The substructuring technique applied for the study presented in 
Section 4 of this report is given in Reference 9, This appendix presents a 
summary of the technique* 

If a structure such as the'one shown in Figure 11 (a) is partitioned 
into substructures, a restrained structure can be defined by assuming just 
those joints at the boundaries of the specified substructure to be 
restrained* Defining the restrained structure in this manner, the various 
substructures shown in Figure 11 (b) become the basic components to be 

used in the analysis of the structural system. The advantage of 
substructuring is that the number of unknown components of joint 
displacements to be evaluated at any one time in the analysis of the total 
system is less than the number of unknowns generated by a direct 
procedure, thus a considerable reduction in CPU time may be realized. 

A substructure can be defined as an element or a group of elements 
that forms a portion of a structural system. The analysis of the 

substructure with its boundaries fixed can be carried out using the direct 
procedure as if it were a complete system in itself. From such an 
analysis, the fixed boundary actions, which are analogous to the fixed end 
actions for a beam element and the displacements of the unrestrained 
joints of the substructure stiffness matrix can also be evaluated. 
Obviously, the joint displacements determined from the analysis will not 
be true true displacements of the unrestrained joints of the substructure 
due to the fact that some of the boundaries of the substructure are 
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the restrained 


artificially constrained in the basic model, i.e., 
structure. Therefore, these displacements must be corrected to account for 
the response of the substructure to the displacement of its artificially 
restrained boundaries, which will be determined from the analysis of the 
total structure. Once the actual displacements of all of the unrestrained 
joints of the total structure have been evaluated, the final end actions 
can be determined for each element of the structure. 


E. 2 Analysis of the Substructure 

For the analysis of the k-th substructure, the static joint 
equilibrium equations can be expressed as follows: 

[ S ] { A } = {P} + {R} E.l 


The various terms in Equation E.l are defined as follows: 

[S] = stiffness matrix, 

{A} = displacement vector, 

{P} = joint load vector, and 
{R} = boundary reaction vector. 

For substructure k. Equation E.l can be partitioned relative to the 
unrestrained (u) and restrained (r) boundary displacements of the 

substructure , as fo Hows : 


nyiiy [M /<yi 

L'Srji^rr'J W-i lip.jj 


+< 


'{ 0 } ■ 

~{\h 


(E.2) 


The true response of a substructure, as part of the total system, is 
the sum of: (i) the response of the substructure, with its boundaries 

restrained, to the external loads; and (ii) the response of the system to 
the actual boundary displacements. The equilibrium equation for the 
response of the substructure due to the displacement of the "restrained" 
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nodes alone can be written as follows: 


Cs tu ']1 [S rr u 


'*ud> 

k'J) 


( 0 }] 
i R rd'J 


(E.3) 


in which {a^}, and {A r( j} are > respectivly, the displacements of the 


unrestrained, and restrained nodes of the substructure due to the 


boundary action vector {R^}, required to maintain equilibrium when the 
bo undries are subjected to the displacement 

Equation E.2 can be expanded as follows: 


[S...J {*.} » {P .J, 

UU U U 

[s ]{A } = {P } + {R }. 
ru u r r 

Solving Equation E.4 gives: 

{A } = [S ] -1 {P } 

U UU u 

substituting which into Equation E.5 leads to: 

{R } - [S ][S ] _1 {P } - {P } 

r ru uu u r 

Similarly, expanding Equation E.3 gives: 

[S ]{A ,} + [S ]{A } = {0}, 
uu ud ur rd 

[S ]{A .} + [S ]{A } = { R } 

ru 1 ud rr ‘ rd rd 

From Equation E.8, 

{A ,} = - [S l'^S ]{A .} 
ud uu ur rd 

substituting which into Equation E.9 leads to: 


(E.4) 

(E.5) 


(E.6) 


(E.7) 


(E.8) 


(E.9) 


(E.10) 


< R r d> - [S s H* rd ) (E.ll) 

in which [S g ] is the substructure stiffness matrix given by: 

[S s ] - [S rr l - [S ru ][s au )' 1 [s iir ] 


E.3 Deflection of Nodes Linking Various Substructures 

The static equilibrium of the nodes linking the various substructures 
k= 1, 2, 3,..., n, and the real foundries can be expressed as follows: 


111 


(E . 13) 


tS a ]{A a } = {P a } + {R a } 


in which the various elements of [S 1 are found as follows: 

a 


(s ). . 

a ij 


<* s A>ij' 


(E. 14) 


where is over k, with k = 1, 2, 3,..., n. Also, {P } is the external 

a 

applied load vector for the nodes linking the various substuctures , and 

{A } is the true deflection vector for the same nodes. {R } is the 
a ' a 

nodal reaction vector* 

Equation E.13 can be written in the following form: 

KiiliM 


l] 

Kdj 

• = - 



J 

,<V~J 


[{P b } J 

1 ■ 


U S 2i^ I ^ S 22^ 

in which { 
generated using the following equation* 
(P rd )i - - + (P) t , 


ff 0 1) 


(E. 15) 


r< j} has been defined in Section E.2. The terms of {P are 


(E. 16) 


where again, 2. is over k, with k = 1, 2, 3,..., n; and {P,} is the 

D 

external load vector corresponding to the actual boundaries of the 
structure. With { . } * {0} Equation E.15 can be expanded to give: 


<A rd ( - ts, , ] -1 {P_ J l, 
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rd 


-1, 


{R b > " I S 21» S nJ {P rd> ' <V 


(E. 17) 
(1.18) 


E.4 Deflection of Unrestrained Nodes 

The true deflection of the unrestrained nodes is obtained as follows: 

k u A ■ «*A + <»ud»k> «- 19 > 

for all k values. Substituting Equation E.10 into Equation E.19 provides: 


<*un>k 


({A } - [S ] _1 [S ]{A ,}).. 

1 u uu ur rd k 


(E.20) 
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APPENDIX F 


FORCE PROGRAM LISTING WITH SAMPLE INPUT/OUTPUT 
FOR SUBSTRUCTURE ANALYSIS OF 3-STORY FRAME 

The information regarding each substructure, the corresponding member 
properties, and the nodal numbering related to the various members are 
read from an input file named SSR.IN6. The input variables for this file 
are SSN, NOM, NOD, NOF, TOR, M, N, MM, NN, L, A, X, where: 

SSN = substructure number, 

NOM =* number of members, 

NOD = total number of nodes times the degree of freedom related to 
each node, 

NOF = number of unrestrained nodes times the degree of freedom related 
to each node, 

TOR = type of member (100 for vertical, and 200 for horizantal), 

M, N =« identification numbers for end 1 of the member, 

MM, NN = identification numbers for end 2 of the member, 

L = member length, 

A = member cross-sectional area, 

X * cross-sectional moment of inertia. 

The identification numbers for the substructures shown in Figure 11 are as 
follows. For the upper substructure , the identification numbers M, N, or 
MM, NN, are: 

•1,3 f o r D , 

4, 6 for E, 

7, 9 for C, and 
10, 12 for F. 
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The identification numbers for the lower substructure are: 


1,3 for B, 

4, 6 for G, 

7, 9 for C, 

10, 12 for F, 

13, 15 for A, and 
16, 18 for H. 

All members have been assumed to be elastic I-sections, with a Young's 
modulus of 29,000 ksi, and the following properties: 

L - 180.0 in. , 

A = 10.3 in. 2 , 

X = 127.4 in. 4 . 

The external nodal loads and the connectivity information about the 
nodes linking the various substructures are read from an input file named 
SSR.IN7. The input variables for this file are: 

LV * the external load vector, 

NODD = the number of boundary nodes in the substructure, 

NS - node connectivity integer between various substructures (0 if 
unconnected, and a non-zero integer if connected). 

The output of the program is directed to a file named SSR.OPT which 
stores the degrees of freedom at the global level, and the corresponding 
deflections. 
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ORIGINAL PAGE 15 
OF POOR QUALITY 


“Ores L b' OT fw :car:t - -- 

: v D gc:,ar5. t i o n o _r a r r a v s a n c v sr: a c 1 a s 

h C. !H 1 — ( C 4 ) , :'i j. CL K vi , O ) , il ( ui , d / i £L \ O , -C .* it ... 

Sparse DOUBLE PRECISION LV ( 18} , wVl < 18) , Ji_RF ( 18) 

3 - a r e c 1 V~ EGER 0 R , 3 , ""CPS. ■ ^ G D , N 0 ~ 4 3 3 N N G M , F , 2 , N s? s . 
3r arse I N~EGE R T I ME 1 , T I ?• EE ? : - L Vf ( £ , Id) 


a 


o n 3. r 8C L/C ---J E< L — H cr. c - d .L U :'v :h : Y l ( £ C ? -£ L ) 


, ij p < 40 ) , 3UU (Id, l 3 ) h 3 R R ( 1 3', 1 3 ) 

IS, 18) 


Shared DOUBLE PRECISION SRU < 13, lfl> , DUF < 1B> , RRF<ia> , S3S C 
Shared DOUBLE PRECISION RRR (18, 18) , 33 (13, IS) , DUF~ t 18) 

Shared DOUBLE PRECISION SS3 <21, £1 ) , RRF3 <21 ) , SLR ( 18, IS) 
Snared DuUdLL P^c.CIdIQN SURr^(13, IS) f S3SS(£, 1 3 5 1 3 ) 

Snared DOUBLE PRECISION SSUR (2, 18, IS) , SSUU <£, 13, 18) , DDL 
Snared DOUBLE PRECISION SSRR<£, IS, 13) , SSRU<£, 18, 18) 
CHPRPC"‘ER*40 F I LE.mm 
E re dec I ar at lore 


< 18, 8) 


Barr i pr 

* ssr, inS-inout file for rnernoer properties data, ssr, in7-iriout file for 

* external nodal loads, ssr. opt-out put file for deflections of the n ocss 

oaen(S, file®* /usr/ul /ravi/ssrS. ir<6 ? . epu— I) 
open ( 1 , f i ie=’ / usr/ u 1 /ravi/ssr5. i n7’ , CPU= 1 ) 
open < £, f i le=^ / usr /u 1 /rav i / ssr5. oot 1 , CPU=1 ) 

CPLL Crrt ic (TIMED 
DO 9 1 = 1, £1 

LV1 <I>=0. 00 
RRF3(I)=0, 00 
DO I® J=i, £1 
3S3 ( I , J ) =0, 00 
10 CONTINUE 
9 CON * I NUE 

End harrier 
Pease 

F I LENM^ 11 /usr/ul /rav i /ssrS. in4 ,J 

Cftww FM (1,4, 1£, 8, SSUU, SSUR, 3SRU, 33RR, FILENM) 

Jsect 

F I LENM= 11 / usr/ u 1 / rav i / ssrS, iri5 ” 

CALL FM <£, 5, IB, 6, SSUU, SSUR, SSRU, S3RR, FILEN*!)^. 

End ocase 

* External nocal loads are read for various su ost ruct ures hy 

* different processors 

Barrier 
oririt*. * OK’ 

NGD=1£ / 

VGP*S 


End :arr 

* E c u a 1 1 c>- n =- £ i s s o 1 v e c u s i r c t n ■= o a r a 1 1 e 1 s u i r o* u 1 1 n e * 3 0 L v 3 
i 'Z- c o r 1 1 r ; Li e 

•L v> •: a 

■X E Hi.) ( ... , H ) ( _ ‘v ( I ) , 1 = 1, -\lCF ) 

Ere carrier 
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•201 Er-.c Drescnec co ORIGINAL PAGE IS 

Prescnsd go 203 1*1, NO? OF POOR QUALITY 

UD £04 w = 1 , r 

£04 SUR ( i , J ) =3311 R < SSN, 1,0 
p 7 ; 3 .End orescnsc 1 co 

-rescnec Do 205 1 = 1 O' 

L/C LLL w = •- , '^Uc 
206 SLLi ( I , J ) =SSULj ( SSN, I , J ) 

205 End oresched do 

Prescned do £08 I = lO 

DO £07 J=1,F 

SRR(I, J)=SSRR(SSN, I, J) 

£07 CONTINUE 

£08 Enc arescnea DQ 


orcsca 1 1 SOLVE (SUU, Lv, DDF, n0;~ ? PV) 

Prescned DO 7 I 1=1, NOF 
DDUr < I , SSN > =DLr ( I > 

PV7(I,SSN>=PV(I> 

71 Ena oresched do 

* Solving equation numoer 3 & 4 using parallel subroutine NPfOlP 

Forcecall PlftTMP (5RU, DDF, RRF, F, NOF, 1) 

Forcecal 1 SGES2 <SUU, SUR, SSS, NOF, F, PV) 

Forcecal 1 MflTMP <SRU, SSS, RRR, F, NOF, F) 

Preached do 117 1 = 1, NOF 
DQ 1 18 J=i, F 
SSSStSSN, I, J) =SSS(I, J) 

118 CONTINUE 

117 End preached do 

Presched DO 81 1 = 1, F 
DO 3£ J=1 , F 

SS(I, J>=SRR<I, J>-RRR<I, J) 

Q£ CONTINUE 
31 End prescned do 

* Forming SRR effective for every substructure and^ placing tne elernen 

* of suostructura 

Barrier 

READ (6,*) NODD 
PRINT*, NODD, ’ NODD GK ? 

READ (6,*) < NS < I), 1=1, NODD) 

P R I NT* , NS < NODD ) , * NS OK 5 
NG= 1 

2 C 1 £ S 1 = 1 = v G D D 

RRF3 ( NO) = RRF3 < NO -OR- ' \G) 

NG=1 

DO £ 1 J= 1 , n.QDD 



C 0 \ ’• \ L ™ 
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vCD= 3 

\QF=S 

Enc oamsr 

" (SSN. SSL i ) 3GT0 : 04 

Barrier 

S3N=3 

NOD = i£ 

NGF=b 
' F=NQD-NGF 
DG 49 1 = 1 , NGF 
LV1 ( I ) = < — 1 • ®) *RRF3 ( I ) 
DO 50 J=l,NOF 
SUU< 1, J) =SS3 ( 1, J) 

3 0 C G N T I N U E 


*+*~f 

CONTI M 

4 

\ 

DO 1 39 
•V ( I ) = 

- Jrf. > t * l J / T \ T — •* 

X , ~ \ — V \ , J. — .f. , 

I*i,NOF 
uV ( I ) +LV1 ( I ) 

1 33 

CONTIN 

UE 


print*, ■ OaK 9 
End barrier 

Forcecai 1 SOLVE (SUU, LV, DUF , NOF, PV) 

For cecal 1 MfiTWS <SSSS, DUF, DUFF, NOF, F, i > 

Barrier 

print*, * OK TILL FERE 9 

* Solving equator 6 for substructure number 1 

DO S9 1=7, 12 

DUFm=DDUF<I-6, l)-DUFF(I-S) 

69 CONTINUE 

End barrier 

Forcecai 1 NATIVE (S3SS, DUF, DUFF, NDF, F , 2) 

Barrier 
print*. 9 OK 9 

* Solving equation S for shstructure nurnoer 2 

DO 79 1=13, 18 

DUF < I ) =DDUF (1-12, 2) -DUF F ( 1-12) 

79 CONTINUE 

l/IhLL CFrt ic ( * I: Y Lu) 
print*, 9 OK 9 
print*. T I MEE, T I : .*E1 

* Writing deflections for common nodes using single processor 

WRITE (2,78) 

WRITE (2,75) 
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DO I*® 1=7, IE 0?, POOR QUAim 

WRITE <£, 77) 1 , DUF ( I ) 

ai R 1 " E ( 2 , * ) 

100 CONTINUE 

* writing deflect ions of noces of suostruct ure II using single processor 
l-I FORItPT (/,SX, » DEFLECTIONS Cr NODES' OF SUBSTRUCTURE 1 1”"/) - 

WRITE < 8,111) 

WRITE <£,75) 

DO 109 1=13, 13 
WRITE (£, 77) I, DUF (I) 

WRITE ( £ , * ) 

109 CONTINUE 

End Banner 

J o i n 

END 

SUBROUTINE FM <SSM, NON, NOD, NGF, SSUU, SSUR, SSRU, SSRR. FI_ENN) 

REAL L, Ki 1 (3, 3) , Kl£ (3, 3) , K2l (3, 3) , K££ (3, 3). 

INTEGER TOR, S3N, F, Z 

DOUBLE PRECISION LV < 18) , LV1 ( 13) 

DOUBLE PRECISION SSUU <£, IS, 13) , SSUR <2, 18, IS) , Sit <13, 13) 

DOUBLE PRECISION SSRU <£, 18, 13) , SSRR (2, 18, 18) 

CHfl ROOTER* A0 FILENM 
OPEN (4, FI LE=F I LENN, CPU=1 ) 

E-29000. 0 

* Sub-structure stiffness matrix formation starts. 

* NOD-Totai degrees of freedom in the substructure, NOF-total nurncer of intern 

* nodes, NOit=number of members, SSN— number of the substructure 

DO 11 1=1, NOD 

DO 12 J=i , NOD ; 

SM<I, J) =0. 8 
18 CONT iNUE 
11 CONTINUE 

* Reads individual member oropert ies, L-Ienath, fl-area of cross-section 

* X— inertia of the section 

* Reads i Y i, N- left si de DOF number ing, XW, NN—ri ght side DOF ^numbering, 

* ™0R is 100 for vertical members, £00 for nors. members 

DC 103 2=1, NON 
READ (4,*) A, X 
Read (4, *) it, n, mit, nn, tor 
crint*, ’ it, N, OK’ 

DO 17 1=1,3 
DO 13 J=i, 3 

J \ - 1 ( j. , J ) =0 
•h «. w ^ 

■<2 1 11, J) =0 
K££ ( I , J ) =0i 

K 11 L, 1 ) = l 2*E*X / ^**3 

\ _ ' * “ , . N } — "G. “V X / h«r ^ v — ) 
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A « * K -t- ^ CL } “K > J. £ CL n *C) / 

K 1 1 (3,3) =A*E*'X/i_ 

*12<1, 1 ) =K 1 1 i i , 1 > * < -l > 

K i CL < 1 , £ ) =K i i ( 1 , 'CL } 

Kl£<l, 3) =K 1 1 (1,3) 

x: 2 (a, d=a:*u,£} 

.•A 1 £ ( £, 2 ) =K 1 1 < a, 2 ) * < — 1 ) 

K 1 E ( cl , ci) = K 11 ( £, 3 / 

•<1£<3, 1 ) =K 1 £ < 1 , 3 ) * < - 1 > 

Kl£ (3, 3> =£*E*X/L 
DO 19 1=1,3 
DO E0 J = 1 , 

K£1 ( I, J> »K1£ ( J, I > 

CONTINUE 

CONTINUE 

K££ ( 1 , i ) =1£*E*X/L**3 
K££ (1,3) =£*E*X/L**£ 

.-<££<£, £ ) =E*A/L 
K££<3, 1 > =6*E*X/L**£ 

K££ < 3, 3) =4*E*X/L 

KK=0 

LL=0 

DO £1 I=MM, NN 

kk=kk+i 

DO ££ J=M, N 
LL=LL+1 

SMCC, J)=S<VUI, J) +K£1 <KK, LL> 

CONTINUE 

LL=0 

CONTINUE 

CONTINUE 

KK=0 

LL=0 

DO £3 I=v*m, NN 
KK=KK+1 
DO £4 J=NM, NN 
!_i_=LL+ 1 

(I , J ) =SM < I , J ) +K££ ( KK , LL ) 
CONTINUE 
LL=0 

CONTINUE 

KK=0 

l_!_=0 

DC £5 I=N, v 
-<K=KK- 1 
DO ££ *VN 

1 

3T (I, J ) =S« < I , J ) +-< 1 £ ( ■-<•< , ) 

CCN“ I >’_■£ 

CDr^IME 


DRK5INAL’ PAGE ES 
OF POOR QUALITY 
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nn 


CRIG3NAE PAGE S3 

23 J=*!, N OF POOR QUALITY 

SN X i , J ) —S?\ X I a J) + K 1 1 XrtK, ■_!_) 

CON" X NLE 

,_ L =0 

CGN~INUE 

GOTO 108 

All < 1, 1) =E*ft/L 

Ki i <2, 2> =12*E*X/L**3 

K 1 1 (2, 3>=6*E*X/L**2 

K 1 i <3, 2) =K i 1 X£, 3) 

K 1 1 (3, 3) =4*E*X/L 
A i 2 X 1 , 1 ) =K 1 i < 1 , 1 > * X - 1 > 

K 1 2 X 2, 2 ) —A 11 <2, 2) *( — 1) 

!*\ _ ^ W ^ """ A L J . X Id , ) 

K12X3, 2)=Kli (2, 3)*< — 1) 

.<12(3, 3 > =£*E*X /' L 
DO 29 I = .1,3 
DO 30 J=l, 3 
K21 (I, J) =K12 X J, I) 

CONTINUE 

CONTINUE 

K 22 < 1 , 1 > = K 1 1 X 1 , 1 > 

K22 (2,2) =K 1 1 (2, 2> 

K22 (2, 3 ) =K 1 2 X 3 , c! ) 

K22X3, 2 ) =K 12X3, 2) 

K22 X 3, 3 ) =K 11 X 3, 3 ) 

KK=® 

LL=0 

DO 31 I=.M, N 
KK=KK+1 ■ 

DO 32 , J=rrtM, NN 
LL=LL+i 

SMXI, J)=SMXI, J) +K12XKK, LL> 

CONTINUE 

LL =0 

CONTINUE 

KK =0 

U _=0 

DO 33 I=N, N 
KK= KK+1 
DO 34 J=N, N 
LL = L . L +1 

3 y i X I, J) =Sr < I, J) — Cl 1 < A' A , > 

_'_=0 

CONTINUE 

KK : =® 


K K •+■ 
38 


\; A 
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=3 

CONTI ML £ 

KK=© 

Li_=0 

DO 37 I=W!vi. 
xk=kk+: 

DO U3 J=N, N 
LL=LL+1 

SM ( I, J)=SM (I, J> +K21 < KK, LL) 

38 CONTINUE 

LL— (S 

37 CONTINUE 

138 CONTINUE 

* Sub structure stiffness matrix formation encs for various 

* sub structures. 

* Natrix for various sue sructures is divided according 

* to number of interior and exterior nodes 

P-NOD-NOF 

DO 1 1 = 1, NICE 

DO 2 J=l, NOF 

SSUU (SSN, I, J) =SM (I, J> 

2 CONTINUE 
1 CONTINUE 

DO 3 1 = 1, F 
DO 4 J = 1 , F 

SSRR (SSiN, I , J) =SM < I+NGF, J+NOF) 

4 CONTINUE 

3 CONTINUE 

DO 5 1=1, NOF 
DO 6 J— 1 , F 

SSUR<SSN, I , J) =SM < I , J+NOF) 

5 CONTINUE 
5 CONTINUE 

DO 7 1 = 1, F 
DO 8 J— 1 , NOF 

SSRU (SSN, I, J)=SM< I+NOF, J) 

S CONTINUE 
7 CONTINUE 
RETURN 
END 

* Subroutine for inversion of matrix (sequetial) 

SUBROUTINE INV<fl,N> 

DOUBLE PRECISION q 
D I TENS I ON S ! 13,13) 
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W| 3 { \i , \ ^ h - ' 

RETURN 

£ND 

* Subroutine for matrix multiplication 

SUBROUT I >E WULT 3, C, *V n> 

DOUBLE PRECISION ft(xa, w(-o) 

DC 3© 1 = 1, N 
C ( I) = 0. 0 
DO 30 L= 1 , N 

30 C < I) =0 < I > + A < 1 , L ) *B < L ) 

RETURN 

* BubrSLs fov matrix muit ioi icM loo (seoiert i*l > 

SUBROUTINE MULTS (A, B, N> _ „ 

DOUBLE PRECISION A ( 13, 13),3<i3, lai , Liid, ' 

DO £0 J=i « N 
DO 30 1 = 1, * 

SUM=0. 0 
DO 40 L= 1 , K 
SUi^SUM+A < 1, L> *8 <L, J> 

40 CONTINUE 
C(I, J)=SUM 
30 CONTINUE 
£0 CONTINUE 
RETURN 
END 

SUBROUTINE MJLT3 IP, B, C, , K, ,N ,U > . 

DOUBLE PRECISION BIS, IB, !S),B(E, .a, *8), Cl. a, -o 

DO £0 J =1 -»N 
DO 30 1=1, M 
SUN=0. 0 
DO 40 L= 1 , K 

3UM=SU!*H"A (II, I , L) *B < 1 1 , L, J > 

40 CONTINUE 

C(I,J)=SUN 
30 CON T INUt 
£0 CONTINUE 
RE T URN 
END 


* Force subroutine 


for matrix multiplication in caral 


-■orcssub Y! ‘”^ ^ 'uti c?« w ? f 


^ -v, \ 4 of ^ 

integer y 2 

3 n a red I -Sc. n — ^ 

Pr i vats I vT~3c.R i * - ? - v __ _ 
i vat a DGUS’-^ ■< ~ - - - - v " ~ ' 

•Z rs c esc. a r' a t c o r • s 

Barrier 

Zr*.c Barr i err*' 

,c -• e sene d D G U * - = - - r ‘ 

EG £03 G 

3...v=0x ? 


lei 
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OF POOt> ottmjty 


SUMeSUF’+ft ( I , <) *B (K, J> 
10® CONTINUE 

C ( I , J ) =sun 
-S 0 CON i I Nub. 

b®0 b T‘ d QrSSCTSC CO 
Barr i er 


End oarrier 
RETURN 
END 


* Force subroutine for solving equation # 7 

Forcesua MftTMS (ft, B, C, M, N, I I ) of NP inert ?*?P 
DOUBLE PRECISION ft < £, IS, 1S),B<1S),C<1S> 
INTEGER m, N, II 
Private INTEGER I,K 
Private DOUBLE PRECISION SUN 
End Declarations 
Barr ier 
End Barrier 
Presched DO 3®0 1 = 1, N 
SUM=0. ® 


DO 1®0 K=1,M 
SUM-SUM+fl < 1 1 , I, K ) *B ( K ) 

100 CONTINUE 
C ( 1 > =SUM 

3®0 End oreschea do 
Barrier 
Ena Barrier 
RETURN 
END 

* force subroutine for solving Cftl Cxi = C bl! in parallel 



C 


Forcesua SOLVE (ft, B, X, N, I PVT) of NP ident NE 

Test program for parallel PLU decomposition using Force versions 
•of "Linoacx" routines. 

T nis program calls a (single stream) routine to read in a matrix 
•anc vector, then SGEFft is called to factor the matrix into a PlU f 
•then SGE3L is called to solve the system. The computation is time 
.anc then a (single stream) routine is called to output the results 


INTEGER N, IPVTU8) 

Shared INTEGER INFO 

DOUBLE PRECISION 0(13,13), B<i8),X<ia> 
Shared INTEGER I TNI, ITVt£, ITW3 
R'-ftu DE'_Tfti, DE-_tqp , RP~E1, Rft~E2 

C G 5 C ci r 5 i 1 O r 3 


Barr i sr 

Chll. C3Vt; ic(I T^l ) 
Zrc; earner 


- orc?ca 1 



33.ii:- Q ( P , 3 - v , 

,V 3 , 0) £;'u 

s a :• re u ar, , * 
sc. i; onre •*, 


-‘V 


: 0 > 


7 rn a “ r i x i s 
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< Yi£ 


44 


C.wJ 


*40 0 

40 1 
4 02 

* ■**■***.< 


,ZUa - ; d£LL Crrt 1C < I7?*£) 

PVT » B, 


-orcscai 1 oGES'L (A, is, ,v, 
c'arr i er 

DC 44 I = i, (v 
X ( I ) =B ( I ) 

call CFrt ic < I TM3 ) 
CALL^PRINTEKA, B, N) 

• ? 


0 ) 


do 25 1 = 1 „ ri 


ELSE 


i ( i. 

J 5 » J = 1 , n ) , 

o(i), x i 

= ( I 

f *Yip __ TT;yn 

) / 50 . 0 

: < 1 

TM 3 - ITM£ 

) /S 0 . 0 

'Ai 

- **E. 0 ) THEN 

— ( 

S. ®*N*N*N/ 

3 . 0 ) * 1 . 


400 ) ,\|P, 

DEL 

<*, 

400 ) NP, 

N, DEL 

<*, 

402 ) 


•END I F 

ir ' ( DEL i A£ , N£. 0 ) THEN 

RftTEE = <N*N)*i. 0E-6 /dp LT A£, 

ELSE ^ '' ~ **’ 401 5 DELTA£ ’ RATES 

u~;I5 <*’ 401 > deltas,' 0.0 

HRxit <*, 402) 

endif 

iCR.MATO Number of processes =* r ■» , w . 

matrix decomoos it ir-, n: » Fl o a , ’ ’ ,a ' cri>< order =* , 14, / 

FORMAT < ’ seconds <’,F5. 3, ’ MFLOPS) » > 

’ rw d/back subst: s rip A , 

FORMAT < * No ai ar , MW ± ■ ’ ’ seconds <*,F5. 3. ’ ,rc:nr«\i« 

End barrier P ** tlme measured.’) ^ ' > 

RETURN 
END 

pc.rc.Bub ssEs*!a,**LDflrc*iP' cr , r*’"****-T************** < ***». 

- "■sv~* o?’ . 

«x-o, using (column = ~ whxCr > solves 

«“* 2 . .«*« ' e :;tr /bacKwar = —■= • 


‘ act or 1 z at a cn o ~ 


: 'V 7 : 


^ovicrss o^ 

wo-<:,ng rca-rix 
J : ~ c sc * arse a i oiens 1 on o f 
tie wording dimension of 
v 1 ® • J ^ v'Ou 

1 © n 3<j. u /out; ;u»,i ~ .. . • .. 

jeo selector's -n *7 7r ~ ?*“ ^ 52£5~„ 

o •; - - j. : so 1 vss "”^ar-sor.=.® t o •• 

^ . Y, iJS r -e c ^. ' -T' " 

rv,n .. _ ^ r, r < Av'.a,-* •» ,. 

-L^L' N “ - - ! V"~ ( i s ) , ;qh "* * ! :: " V5 

•;. -3 -- ’— ■ v ' «. O r- 1 » I , 4 ; 

? :r--cs:pv' ;;; „ 


- o l. s 


124 



ORIGINAL PAGE IS 

Sparse DG,.3_E PRECISION ?S>| OF POOR QUALITY 

Ere ceclaratioris 


C first: solve: i_*v = 3 

CO £0 K = X , Ns-X 

Barr i er 

L = I PVT (K) 

TEM = B(L) 

B < i_) = B < K > 

B<K) = TEM 
End barrier 

Prescnea do ££ I = K+l, N 
B ( I ) = B <. I ) + TEM * 3!2,K) 
££ End orescned do 

20 continue 


C then solve: U*X = Y 

do 42 K = N, 1, -X 
Barrier 

BOO = BOO /fi(K, K) 

TEM = -B(K) 

End barrier 

Prescned do AS 1=1, K-i 
B < I ) = B < I > + TEM * fl(I,K) 

AS End oresehed do 

A0 continue 

Barrier 
End barrier 
RETURN 
END 

********-**##***•#•*****#*******#****#****************.*******#***#*#***•**•*•* 
Forces ub 3GEFA(A, LDA, N, IPVT, INFO) of NP ident ME 
INTEGER LDA, N, IPVT (LDA), INFO 
DOUBLE PRECISION A (LDA, LDA) 

C 

Private DOUBLE PRECISION ABSMAX 
Private DOUBLE PRECISION T 
Private INTEGER XK, I, MAXI, J, KP1 
Shared INTEGER IA:_^ 

Shared DOUBLE PRECISION PIV, ALLMAX 
Shared logical I LOCK 
Ena declarations 


Barr 1 er 


1 nF G ■= a 

■~u w M X = -i> m 


;-t s m a •" x 


1000 « - 1, V-X 

(InFG. N£. 8) GG T 0 2000 


; e s et . o cj . -<i a x •>. ,r. a 
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FBStfAX = 0. 0 

C. • Find the o i vot row 

: ‘rsscned DC 10® 1 = KK. N 

T = AB3( A ( 2 , KK) ) 

IF (ABSMAX . LT. T) THEN 
HOAX I » I 
ABSMAX = T 
ENDIF 

1 0® End Presched DC 

Cm Locate the Sharec Information if riscassary 

Critical I LOCK 

IF ( ALLMAX . LT. AB3MAX) THEN 
I ALL. » TAX I 
ALLMAX = AB3MAX 
ENDIF 

End critical 
Barrier 

I h ( ALLMAX . EQ. Q}» 0) THEN 
INFO = KK 
I PVT <KK> =KK 

ELSE 

I PVT ( KK ) = I ALL 
END IF 

Ena Barrier 

C ..If the matrix is singular then oass the information 


IF (INFO . EQ. 

KK) 

GOTO 1000 

MAX I = I ALL 



IF (MAXI . NE. 

KK) 

THEN 


C. ...... . Swao rows if necessary 

Preserved DO 11® J » KK, N 
TEMP = A (MAX I, J) 

A (MAX I , j ) = P(KK, J) 

A (KK, J ) = T~*P 
End -resened DG 
END I r 

D. . • . . « . Self schedule row reductions 

KPl = KK - 1 
Barr i sr 

I V = ~ }. , 0 / d ( •< r< , K X ) 


= I V * A ( 
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pc.k-o = 

di - C! l&* o ■“ K 4 I , \\ 

f-i < .1 , J > — A<I, J) + <KK, J) 

- — ® C 0 i\ T I N U E 

ZnG Seifsc'ieo DG 

Barr x er 

ALUNPX=0 
End barrier' 

i®0® CONTINUE 

2 000 RETURN 

END 

-orcesuo 3GES2 < A, B, X, N, in, I-'vT> of NP icem ,r 
INTEGER IPVi < IS) , JGB, :V 

DGuBuE PRECISION A ( 13, IS) , B < IS, IS) , X < 18, 13) 

Private integer I , K, L, J 

Snared DCUBuE PRECISION Tr>1 

End deciarat ions 

DG 30 J = 1 , N 

♦ first solve i_*V=B 

DO £0 K— 1 , N— 1 
Barr x er 
L— 1 PVT ( K ) 

TM=B < L, J) 

B(L, J) =B (K, J) 

3 (K, J) *TM 
End oamer 

Preschsd go £2 I=K+1,N 
3 < I , J ) =B (I , J ) +Ti*l*A < I , K > 

E£ End preschea do 

£0 continue 

* then solve U*X=y 

DO 40 K=N, 1,-1 
Barrier 

3(K, J) =B ( K, J) K) 

Tjvj=-B <K, J ) 

Enc barrier 

Preschea ao 4£ 1 = 1 , K-l 
B < I , J) =8 < I , J) +7M* A ( I , K ) 

AS Enc oresched do 

A 0 c o yr,Ji z i u © 
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ynaresP r '&09BLS -'RECiSICN ”N 

End c&d'CaeaCions 

DO 30 

* EIRST SOLV7S =r*y»B 
DO 20 N-l 

Barrier ■ ' 

i_= : PVT CK , II) 

TN=B < 1 1 , L, J) 

Bill, L, J)=B<II, K, J) 

3 < I 2 , K, J > =TN 
End barrier 

Prescned DO ££ I=K+1,N 
3< II, I, J> =B< : 2, I, J) +TN*fl< : I , I, K> 
2£ End Prescnea do 
£0 CONTINUE 
♦MOW SOLVE U*X = Y 

DO 40 K=N, 1, -i 
Barr ier 

B(II,K,J)=B(II,K, J) / ft < 1 1 , K , K > 
TM=-B < 1 1 , K, J) 

End barrier 
Prescned do 4£ 1=1, K- 1 
B ( 1 1 , I , J) =B (II, I, J) +TN*fl(II, I , K ) 
42 End Presched do 

40 CONTINUE 

DO 50 1=1, M 
C<I, J)=B<II, I, J) 

50 CONTINUE 

30 CONTINUE 

RETURN 
END 
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ssr'. xri6 

1 2 , £ , 4 , 1 

130. 0, 10. 3, i£7, 0 

a. , , 4 . & . £. 0 0 

130. 0, 10. 3, 127. 0 
7, 3, 10, 12, 200 
130. 0, 10. 3, 127. 0 
7,9, 1,3, 100 
130. 0, 10. 3, 127. 0 
10, 12, 4, 2, 100 

.L d 

x, to * 4, ■ — i , to , 0, 0, 0, 0, 0, 0 

13,2,5,2 

130. 0, 10. 3, 127. 0 
1 , 3, 4, 2, 200 

180. 0, 10. 3, 127, 0 
1, 3, 7, 9, 100 

180. 0, 10. 3, 127. 0 
4, 2, 10, 12, 100 
180. 0, 10. 3, 127. 0 

13, 15, 1, 3, 100 
180. 0, 10. 3, 127. 0 
12, 18, 4, 2, 100 

12 

1, 2, 3, 4, 5, 2, 7, 8, 9, 10, 11, 12 
7 
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JOB = 


>sr. in/ 


3sr« 


ir»7 


■-Z. C' a ? c‘ * “ 

J . d , 0.0 

is. a, a. 

3 C$2 


3 , 0 , 0 , 0 . : 2 ? ®* ®? ® 

0, 0. 0, a. a? 0 - ®» 0* 0 
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5 Sr'- 0 3t 

"’"IONS OF CODE'S CO'" 3 

ORIGINAL PAGE IS 

of poor quality: 

*nQN to SUBSTRUCTURES 

i. CF freEDGM 

DEFLECT :DN 

1 

1 1 . 72159 

a 

0. 06337 

3 

- 0. 02147 

jLl 

11.71556 

3 

- 0. 06337 

b 

- 0. 02145 


DEFLECTIONS of nodes 

OF SUBSTRUCTURE I 

DEG. OF FREEDOM 

DEFLECTION 

7 

15. 6S023 

a 

0. 07636 

9 

-0. 01006 

10 

i5. 67420 

1 1 

-0. 07 696 

IE 

-0. 01003 

DEFLECTIONS OF NODES 

OF SUBSTRUCTURE II 

DEG. FREEDOM 

DEFLECTION 

3 

4. 33344 

14 

®. 04478 

1 - 

-0, 02506 

l S 

4. 89343 

‘ 7 

—0. 044,- 3 

i a 

—0i. 025©*- 
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APPENDIX G 


TITLES OF ABSTRACTS FOR 29TH AIAA/SDM CONFERENCE 


The following abstracts are submitted for the 29th AIAA/SDM 
Conference to be held in Williamsburg, Virginia, April 18-20, 1988: 


1. Prasad, V., Razzaq, Z. , and Storaasli, 0.0., " Two-Dimensional 

Concurrent Finite Element Analysis of Rectangular Panel with Hole." 

2. Bhati, R. , Razzaq, Z. , and Storaasli, 0.0., "Concurrent Substructure 
Analysis of Plane Frames." 
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